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ABSTRACT 


A  multichannel  random  accessed  communication  system  was 
considered.  Every  time  slotted  channel  is  shared  by  a  large 
number  of  transmitters.  To  solve  a  collision  which  occurs  at 
a  certain  time  slot,  several  collision  resolution  proto¬ 
cols  were  tested,  analyzed,  and  simulated.  The  quantities 
to  be  investigated  in  this  thesis  include  (the  average 
number  of  slots  required  to  resolve  a  collision  involving  n 
users),  linear  bounding  of  and  the  maximal  achievable 
throughput . 
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I.  INTRODUCTION 


The  purpose  of  this  thesis  is  to  propose  and  analyze  a 
few  collision  resolution  protocols  suitable  for  multichannel 
environment.  To  have  some  familiarity  with  collision 
resolution  algorithms,  the  computer  communication  networks 
and  the  progression  of  the  algorithms  should  be  explained. 


A.  BACKGROUND 

1 .  Computer  Networks 

In  the  past  several  decades,  we  have  witnessed  the 
birth  and  unprecedented  growth  of  the  computer  industry.  In 
this  spectacular  progression  of  the  computer  industry, 
computers  have  become  smaller,  cheaper  and  more  numerous. 

The  increased  use  of  computers  in  the  communications 
field  has  greatly  influenced  the  design  and  architecture  of 
computer  systems  and  networks.  In  the  past,  computing  was 
generally  performed  at  a  centralized  facility,  using  a 
mainframe  computer.  Present  day  systems  increasingly  rely 
on  distributed  processing  using  smaller  remotely  located 
computers  connected  by  a  network.  These  systems  are  called 
"computer  networks".  A  computer  network  can  be  defined  as 
an  interconnected  group  of  independent  computer  systems 
which  communicate  with  one  another.  These  computer  systems 
share  resources  such  as  programs,  data,  hardware,  or 
software.  Obviously,  the  goals  of  computer  networks  are 
resource  sharing,  load  balancing,  high  reliability,  and  the 
reduction  of  communication  costs.  Finally,  it  can  provide  a 
powerful  communication  medium  among  widely  separated  people. 


Figure  1.1  shows  the  conceptual  model  of  a  computer 
network.  Two  basic  components  of  the  communication  subnet 
are  switching  elements  and  transmission  lines.  The 
communication  subnet  is  responsible  for  transmitting 
messages  reliably  and  rapidly  between  all  participating 
network  sites. 


Figure  1.1  Computer  Network 


2 .  Computer  Communication  Networks 

A  computer  network  consists  of  communication 
devices,  the  host  processors,  the  transmission  lines  and  a 
set  of  rules.  These  rules  can  be  implemented  in  either 
hardware  or  software  to  insure  the  orderly  flow  of  traffic 
in  the  network.  As  computer  networks  are  expanded, 
communication  between  computers  has  become  more  complicated 
and  more  important.  For  data  communication  in  a  computer 
communication  network,  several  switching  techniques  can  be 
used.  These  are  circuit  switching,  message  switching,  and 
packet  switching. 


Circuit  switching  is  very  similar  to  telephone 
conversation  and  once  a  route  has  been  set  up,  it  will  be 
used  for  the  entire  conversation.  Message  switching  is 
another  switching  strategy.  A  transaction  between  two  users 
in  the  network  can  be  visualized  as  a  conversation  and  each 
message  is  transmitted  independently  in  the  manner  of 
store-and-forward.  In  packet  switching,  each  message  breaks 
into  several  fixed- length  data  units  called  packets.  Main 
advantages  of  packet  switching  include  reduced  delay  and 
improved  throughput.  It  is  particularly  suitable  in  handling 
interactive  traffics  such  as  data  generated  by  computers. 
Therefore,  computer  communication  is  usually  done  on  the 
basis  of  packet  switching. 


3 .  Problems  Involved  in  Computer  Communication  Network 

There  are  several  problems  that  should  be  considered 
in  designing  a  computer  communication  network.  These 
include  network  topology,  routing  algorithm,  congestion 
contro.l,  capacity  allocation,  concentrator  siting  and 
network  access  techniques. 

In  this  thesis,  network  access  techniques  will  be 
studied.  It  is  desirable,  if  interference  can  be 
suppressed,  to  transmit  and  receive  a  number  of  signals  over 
a  single  communication  path.  By  sharing  the  same 
communications  path  with  other  signals  adopting  such  a 
"multiple  access"  technique  may  effectively  reduce  the  cost 
per  bit  of  data  transmitted. 

4 .  Multiple  Access  Techniques 

One  classical  solution  to  the  multiple  access 
problem  is  multiplexing,  e.g.,  time  division  multiple 
access  (TDMA)  or  frequency  division  multiple  access  (FDMA) . 
In  FDMA,  all  users  share  the  frequency  spectrum  of  a 
communications  path  and  transmit  simultaneously.  Note  that 


each  user  is  allocated  a  unique  frequency  band.  If  the  full 
bandwidth  W  of  the  communications  path  is  divided  into  N 
users  (or  N  channels),  then  each  user  can  transmit  only  at 
speeds  less  than  the  frequency  slot  of  W/N.  The  limitation 
is  due  to  the  need  for  small  guard  bands  between  adjacent 
channels  to  prevent  any  sideband  signal  from  overlapping.  In 
TDMA,  all  users  occupy  the  same  RF  bandwidth,  but  transmit 
sequentially  in  time. 

The  next  solution  of  the  multiple  access  problem  is 
random  access  technique,  which  permits  any  user  to  seize  the 
entire  communications  resource  when  he  has  information  to 
transmit.  Random  multiple  access  is  one  of  the  most 
efficient  tools  for  channel  sharing  by  a  number  of  users. 
Random  access  technique  is  expedient  when  there  are  a  large 
number  of  stations  with  bursty  input  traffic. 

5 .  Random  Access  Techniques 


In  1970,  Abramson  [Ref.  14]  devised  a  pure  ALOHA 
system.  In  this  system,  a  single  channel  is  shared  by  all 
users  and  a  user  transmits  data  as  soon  as  he  is  ready  to  do 
so.  It  is  assumed  that  the  data  is  transmitted  in  packet 
form  and  all  transmitted  packets  have  identical  length. 
Each  packet  contains  the  information  to  be  transmitted, 
along  with  the  source  and  destination  addresses  and  anv 
other  overhead  information  (error  correcting  and  error 
detecting  bits).  It  is  supposed  that  the  transmitter  only 
sends  a  packet  of  some  fixed  duration,  say  T  seconds.  During 
this  T-second  transmission,  if  only  one  transmitter  sends  a 
packet,  then  the  transmission  is  successful.  Since  there  is 
no  synchronization  among  users,  transmission  times  of  two  or 
more  users  may  overlap.  This  situation  is  referred  to  as  a 
collision.  However,  due  to  the  perfect  feedback  property  of 
packet  broadcasting,  the  sender  of  a  packet  can  always  find 
out  whether  or  not  his  packet  was  successfully  transmitted 


by  Just  listening  to  the  channel.  A  collision  results  in 
complete  loss  of  the  information  carried  by  the  collided 
packets,  thus  retransmission  of  such  packets  are  necessary. 
To  avoid  a  repetition  of  the  same  collision,  pure  ALOHA 
specifies  that,  after  a  collision,  each  transmitter  involved 
randomly  selects  a  waiting  time  before  it  again  retransmits 
its  packet.  Throughput  of  the  system  can  be  measured  as  the 
average  number  of  successful  transmissions  during  a  message 
transmission  time  (usually  normalized  to  1  for  convenience), 
and  is  given  by: 

S  =  G  Po 

where  S  is  the  throughput,  G  is  the  channel  load  measured  as 
the  number  of  attempted  transmissions  per  unit  of  normalized 
time,  and  Pg  is  the  probability  of  no  collision.  Assuming 
that  the  channel  load  G  follows  Poisson  distribution,  P„  for 
pure  ALOHA  is  given  by: 


therefore: 

S  =  Ge'^^ 

The  maximum  throughput  one  can  have  is  0.184, 
corresponding  to  G=0.5. 

In  1972,  Roberts  [Ref.  15]  modified  Abramson’s 
algorithm.  This  algorithm  is  now  known  as  the  slotted-ALOHA 
random  access  algorithm  which  doubles  the  capacity  of  the 
pure  ALOHA  system.  This  scheme  operates  exactly  the  same  as 
the  pure  ALOHA  does,  except  that  a  new  packet  transmission 
must  begin  at  the  next  slot  boundary.  Time  is  divided  into 
slots  with  length  equal  to  the  packet  transmission  time.  If 
a  packet  is  generated  in  the  middle  of  a  slot,  then  it  must 
wait  until  the  next  slot  boundary  before  it  is  transmitted. 
A  collided  packet  must  be  retransmitted  at  some  random 


.TO  < 


future  time  slot.  The  probability  of  no  collision  in  this 
case  is  given  by: 

and  S  =  Ge'®  gives  a  maximum  throughput  of  0.368  for  G=l. 

Small  increases  in  the  channel  load  can  drastically 
reduce  its  performance  because  the  number  of  collisions 
increases  exponentially  with  higher  values  of  G.  It  should 
be  noted  that  system  performance  can  be  evaluated  by  three 
properties:  throughput,  delay,  and  stability.  ALOHA  channels 
are  fundamentally  unstable,  but  there  exist  a  number  of 
simple  control  procedures  which  stabilize  these  channels. 
Most  of  these  control  schemes  estimate  the  number  of  busy 
terminals  in  the  system.  The  "Collision  Resolution 
Approach"  opens  the  possibility  of  an  exact  analysis  of  the 
stable  steady-state  behavior  of  a  random-access  system. 

B.  HISTORICAL  REVIEW  OF  COLLISION  RESOLUTION  ALGORITHMS 

Over  the  years,  the  research  on  collision  resolution 
algorithms  has  always  focused  on  the  bursty  arrivals  of 
messages  and  the  interference  between  transmitters,  but  has 
generally  ignored  the  noise.  It  has  been  always  assumed 
that  a  message  can  be  successfully  communicated  in  the 
absence  of  collision.  A  Collision  Resolution  Algorithm  (CRA) 
can  be  defined  as  an  algorithm  which  organizes  the 
retransmission  of  colliding  packets  in  such  a  way  that  every 
packet  is  eventually  transmitted  successfully  with  finite 
delay.  Collision  resolution  interval  (CRI)  is  referred  to  as 
the  period  beginning  with  the  slot  containing  the  original 
collision  (if  any)  and  ending  with  the  slot  in  which  the 
original  collision  is  resolved  (or  ending  with  the  first 
slot  when  that  slot  is  collision- free ) .  General  assumptions 
for  the  CRA  are  as  follows. 


The  forward  channel  to  the  receiver  is  a  single, 
errorless,  time-slotted,  collision-type  channel. 

The  transmitter  can  send  a  data  in  a  packet  form  whose 
duration  is  same  as  one  slot. 

A  "collision"  between  two  or  more  packets  is  always 
detected,  at  least,  by  the  participating  receivers. 

A  noiseless  feedback  broadcast  concerning  the  status  of 
a  slot  is  always  available  to  all  the  participating 
users.  The  status  of  a  slot  could  be  any  one  of  the 
following:  (a)  that  slot  was  empty  or  («)  that  slot 


two  or  more  packets. 

Signal  propagation  delays  can  be  ignored. 


.  Capetanakis  Collision  Resolution  Algorithm 


The  recent  discovery  by  Capetanakis  [Ref.  6]  of  a 
collision-resolution  algorithm  was  a  surprising  development 
in  the  evolution  of  random-access  techniques.  This  tree 
algorithm  shall  be  referred  to  as  the  Capetanakis  Collision 
Resolution  Algorithm  (CCRA).  The  CCRA  can  be  stated  as 
follows:  after  a  collision,  all  transmitters  involved  flip 
a  binary  fair  coin,  those  flipping  0  retransmit  in  the  very 
next  slot,  those  flipping  1  retransmit  in  the  next  slot 
after  the  collision  (if  any)  among  those  flipping  0  has  been 
resolved.  No  new  packets  are  allowed  to  transmit  until  after 
the  initial  collision  has  been  resolved  [Ref.  8:p.  77]. 

The  following  example  should  both  clarify  the 
algorithm  and  illustrate  its  main  features.  Suppose  that 
the  initial  collision  is  among  three  users,  as  shown  in 
Figure  1.2.  For  convenience,  we  refer  to  these  users  as  A, 
B  and  C.  After  the  collision  in  slot  1,  all  three  of  these 
users  flip  their  coins;  we  suppose  that  C  flips  0  while  A 
and  B  flip  1.  Thus  C  sends  in  slot  2,  while  A  and  B  transmit 
in  slot  3,  and  a  collision  occurs.  After  a  collision  in  slot 
3,  we  suppose  that  A  and  B  both  flip  1.  Thus,  slot  4  is 
empty,  so  A  and  B  again  recognize  that  they  should 
retransmit  in  slot  5.  After  a  collision  in  slot  5.  we 


suppose  that  A  and  B  both  flip  0.  Hence  they  both  retransmit 
in  slot  6.  After  the  collision  in  slot  6,  we  suppose  that  A 
flips  0  and  B  flips  1.  Thus,  A  successfully  transmits  in 
slot  7  and  B  successfully  transmits  in  slot  8.  All  three 
users  in  the  original  collision  have  now  transmitted 
successfully,  but  the  collision  is  not  yet  resolved.  The 
reason  for  this  is  that  no  one  can  be  sure  that  there  was 
not  another  user  who  transmitted  in  slot  1,  then  flipped  1 
in  collision  slot  3  and  5  and  who  thus  is  now  waiting  to 
retransmit  in  slot  9.  It  is  not  until  slot  9  proves  to  be 
empty  that  the  original  collision  is  finally  resolved. 

Figure  1.3  shows  the  tree  diagram  which  illustrates 
the  same  situation  as  Figure  1.2.  The  number  inside  the 
node  indicates  the  feedback  information  for  that  time  slot 
(0  =  empty  slot,  1  =  single  packet,  >2  =  collision),  and  the 
binary  numbers  on  the  branches  coming  from  a  node  indicate 
the  path  followed  by  these  users  that  flipped  that  binary 
number  after  the  collision  at  that  node.  The  number  above 
each  node  indicates  the  time  slot.  At  slot  9,  all  users  now 
simultaneously  learn  that  all  the  packets  in  the  original 
collision  have  been  successfully  transmitted  because  a 
binary  rooted  tree  has  two  more  terminal  nodes  than  it  has 
intermediate  nodes  excluding  the  root  node.  This  is  the 
property  of  CCRA.  Additionally,  all  transmitters  must  know 
when  the  original  collision  (if  any)  has  been  resolved  as 
this  determines  when  n^w  packets  may  be  sent. 

In  the  implementation  of  CCRA,  the  feedback 
information  is  used  only  to  determine  whether  the 
corresponding  slot  has  a  collision  or  is  collision- free . 
Therefore,  it  is  not  important  whether  feedback  information 
is  empty  or  successful  (only  one  packet),  as  both  cases  mean 
that  the  slot  is  collision-free. 
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Figure  1.2  Collision  Resolution  Interval  for  CCRA 


Figure  1.3  Tree  Diagram  for  CCRA 

The  main  interest  is  the  length  of  the  CRI ,  i.e.,  of 
the  number  of  slots  in  the  CRI,  because  system  throughput  is 
determined  by  dividing  CRI  into  the  number  of  packets 
involved  in  the  original  collision.  Capetanakis  proves  that 
his  scheme  achieves  throughputs  above  the  1/e  "barrier"  for 
slotted-ALOHA.  The  throughput  for  the  Poisson  source  model 
and  an  infinite  number  of  independent  identical 
sources  (users)  is  0.346.  He  also  proposed  a  dynamic  CRA 
which  leads  to  a  maximum  throughput  of  0.430.  A  similar 
idea  was  proposed  by  J.  F.  Hayes  [Ref.  7]. 
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2 .  Massey’s  Algorithm 


Observing  that  if  a  collision  slot  is  followed  by  an 
empty  slot,  one  slot  can  be  saved  by  repeating  the  random 
retransmission  (realized  by  probability  0.5)  before  a 
predictably  certain  collision  is  allowed.  Massey  [Ref.  8] 
improved  the  Capetanakis  algorithm,  called  Modified 
Capetanakis  Collision  Resolution  Algorithm  (MCCRA) .  MCCRA 
is  like  the  CCRA  algorithm,  except  when  the  feedback 
indicates  a  slot  in  which  a  set  of  users  who  flipped  0 
should  retransmit  is  empty.  Then  each  user  involved  in  the 
most  recent  collision  flips  a  binary  fair  coin.  Those 
flipping  0  retransmit  in  the  very  next  slot,  those  flipping 
1  retransmit  in  the  next  slot  after  the  collision  (if  any) 
among  those  flipping  0  is  resolved.  The  following  example 
shows  how  MCCRA  works. 

Figure  1.4  gives  a  binary  tree  for  a  CRI  containing 
two  packets  in  the  original  collision,  and  for  which  both 
users  flipped  1  on  the  first  two  tosses  of  their  binary 
coin,  but  one  flipped  a  0  and  the  other  a  1  on  their  third 
toss.  The  nodes  labelled  "skip”  and  having  no  slot  number 
written  above  them  correspond  to  points  where  the  feedback 
indicates  that  certain  users  should  immediately  flip  their 
binary  coins  to  thwart  a  certain  collision.  In  this  example, 
we  need  two  fewer  slots  than  CCRA  because  the  nodes  labelled 
"skip”  would  become  collision  slots  [Ref.  8:p.  80]. 

MCCRA,  unlike  CCRA,  requires  the  feedback 
information  to  distinguish  between  empty  slots  and  the  slots 
with  one  packet.  Massey's  modified  algorithm  (MCCRA)  induces 
a  throughput  equal  to  0.375  in  the  static  case,  and  a 
throughput  equal  to  0.462  for  the  dynamic  case.  However, 
MCCRA  algorithms  are  very  sensitive  to  channel  errors.  If 
channel  errors  occur,  the  CRI  never  terminates  and  no 
packets  are  transmitted  after  some  point. 


Figure  1.4  Tree  Diagram  for  MCCRA 


3 .  GPK  Algorithm 

In  1982,  Georgiadis  and  Papantoni-Kazakos  [Ref.  9] 
proposed  a  further  improved  CRA  which  is  stable  at  higher 
input  rate.  Their  algorithm  involves  the  knowledge  of 
multiplicity,  i.e.,  the  number  of  users  involved  in  a 
transmission.  At  first,  they  assumed  that  the  number  of 
packets  involved  in  a  collision  can  always  be  detected 
correctly  by  a  device  called  "energy  detector".  This 
protocol  is  called  a  Collision  Resolution  with  Additional 
Information  (CRAI).  CRAI  has  a  first-corae-first-serve 
feature  and  can  be  described  as  follows. 

Suppose  n  users  are  known  to  be  involved  in  a 
collision.  Then  each  of  these  users  independently  flips  a 
bias  coin  with  P(H)  =  0^  and  P(T)  =  l-Ofj  to  break  into  two 
subgroups.  Collision  (if  any)  in  the  second  group  will  not 
be  resolved  until  all  users  in  the  first  group  have 
successfully  transmitted.  The  optimal  value  of  can  be 
determined  by  classic  dynamic  programming. 

Again,  they  proposed  a  modified  version  of 
CRAI  (MCRAI)  for  the  case  in  which  the  energy  detector  is 


capable  of  telling  the  multiplicity  of  a  transmission  up  to 
a  certain  limit.  This  MCRAI  maintains  the  FCFS 
characteristics  of  CRAI  and  has  certain  robust  properties. 
It  means  that  MCRAI  is  robust  in  the  presence  of  certain 
channel  errors. 

Georgiadis  and  Papantoni-Kazakos  observed  that  CRAI 
stability  region  is  obtained  for  an  input  rate  below  0.5324. 
Stability  region  of  MCRAI,  for  energy  detectors  of  capacity 
7,  is  very  much  the  same  as  the  stability  region  of  CRAI. 

Almost  the  same  type  of  CRA  was  developed 
independently  in  the  USSR  by  Tsybakov  and  Mikhailov  [Ref. 
11] .  A  slightly  different  conflict  resolution  algorithm  was 
introduced  by  Gallager  [Ref.  10]  for  the  same  model  as  in 
Capetanakis'  and  Massey's  models.  His  algorithm  also  has 
stability  and  FCFS  characteristics.  Also,  the  algorithm 
realizes  its  maximum  throughput  at  0.4872.  Gallager' s 
protocol  was  later  improved  by  Humblet  [Ref.  12],  with  an 
increased  throughput  0.48775. 

C.  SCOPE  OF  RESEARCH 

As  seen  in  previous  sections,  recent  research  on  random 
access  techniques  has  been  centering  around  efficient 
collision  resolution  algorithms.  So  far,  most  of  the  work 
reported  in  the  literature  only  consider  single  random 
accessed  channel.  Collision  resolution  algorithms  applicable 
to  a  multichannel  system  have  rarely  been  studied.  The 
purpose  of  this  thesis  is  to  propose  and  analyze  a  few 
collision  resolution  protocols  suitable  for  multichannel 
environment.  For  the  multichannel  CRA's,  one  can  group  them 
into  the  following  categories. 

Category  1  -  The  collision  in  each  channel  is  resolved 
independently.  New  transmissions  should  not  be  attempted 


until  collision  in  all  channels  are  cleared.  At  the 
beginning  of  a  CRI,  a  station  (user)  with  a  packet  to  be 
transmitted  randomly  picks  one  channel  and  transmits.  In 
selecting  users  for  transmissions  at  the  beginning  of  a  CRI, 
Gallager's  idea  [Ref.  10]  can  be  used.  All  ready  users  may 
also  be  allowed  to  transmit  at  the  beginning  of  a  CRI. 

Category  2  -  Each  channel  runs  independently.  At  the 
beginning  of  the  CRI  of  a  certain  channel,  we  can  either 
allow  all  ready  users  to  transmit  or  select  users  for 
transmission  using  Gallager's  idea. 

Category  3  -  Collisions  in  all  channels  are  resolved 
jointly  by  efficient  algorithms.  Otherwise  it  remains 
identical  to  the  first  category. 

In  this  thesis,  category  3  is  emphasized,  and  some 
collision  resolution  protocols  are  proposed.  It  is  further 
assumed  that  channel  status  information  on  the  feedback 
channel  is  either  ternary  or  m-ary.  A  ternary  feedback 
specifies  whether  the  channel  is  empty,  successful,  or 
colliding.  An  m-ary  feedback  adds  the  information  about  the 
number  of  users  participating  in  a  transmission.  Under 
category  3,  the  following  three  collision  resolution 
algorithms  are  studied.  Detailed  descriptions  will  be  given 
in  each  corresponding  chapter. 

1 .  Static  M-ary  Splitting 

At  time  t,  each  of  the  packets  arriving  during  in 
"^0’  '■0*  transmitted  through  a  randomly  picked 
channel.  Packets  arriving  in  "tQ-^  A  .  tg-f  2A"  will  not  be 
transmitted  until  any  collisions  of  the  previously 
transmitted  packets  are  completely  resolved.  Energy 
detectors  are  not  available.  When  collision  is  detected  at 
a  certain  time  slot,  all  packets  involved  in  collision 
within  this  time  slot  are  divided  into  m  groups.  Each  packet 
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is  randomly  placed  in  one  of  m  groups  by  flipping  an  m- faced 
dice.  The  packets  in  the  first  group  are  retransmited 
immediately  through  one  of  the  randomly  picked  available 
channels.  The  second  group  is  allowed  to  transmit  after  all 
of  the  packets  in  the  first  group  are  successfully 
transmitted,  and  so  on. 

2 .  Dynamic  M-ary  Splitting 

An  infinite  energy  detector  is  made  available  in  the 
static  m-ary  splitting  protocol.  Therefore,  the  number  of 
packets  involved  in  a  collision  is  known  and  we  can 
accordingly  determine  the  value  of  m  (number  of  groups) 
which  minimizes  the  CRT.  The  rest  of  the  protocol  remains 
the  same  as  in  static  m-ary  splitting. 

3 .  Binary  Splitting  with  Energy  Detector 

Assume  the  availability  of  an  infinite  energy 
detector.  After  a  collision,  each  user  involved  in  the 
collision  randomly  flips  a  biased  coin  to  either  join  the 
first  or  the  second  group.  Let  denote  the  probability 
that  a  user  will  join  the  first  group  when  the  collision 
involves  "n"  users.  Again,  we  shall  optimize  To 

transmit  a  packet,  one  of  the  channels  again  will  be  picked 
randomly.  If  a  user  decides  not  to  retransmit  immediately, 
then  its  packet  will  be  held  until  any  collision  involved  in 
the  previous  slot  is  resolved. 

D.  OUTLINE  OF  THIS  THESIS 

In  Chapter  II,  static  m-ary  splitting  is  analyzed  and 
compared  with  simulation.  The  main  results  include: 


a)  The  average  number  of  slots,  L^^,  required  to  resolve  a 
collision  involving  "n"  users  is  obtained.  Through 
numerical  examples,  it  is  observed  that  binary 
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splitting  achieves  the  minimum  L^j.  The  values  of 
are  verified  by  computer  simulation. 

b)  It  is  proved  that  is  linearly  bounded.  That  is, 

there  exist  n^, ,  ‘  i^un  such  that  -C  < 

^iSuno'^  n  >  n^  . 

c)  Based  on  the  linear  bounding  of  L_  the  maximum 
achievable  throughput  is  obtained. 

Dynamic  m-ary  splitting  is  considered  in  Chapter  III. 
Results  similar  to  those  mentioned  above  are  obtained. 
Dynamic  m-ary  splitting  performs  better  than  the  best  static 
m-ary  splitting,  i.e.,  the  static  binary  splitting  studied 
in  Chapter  II. 

In  Chapter  IV,  the  binary  splitting  with  energy  detector 
protocol  is  analyzed.  It  is  again  observed  that  the  dynamic 
m-ary  splitting  outperforms  this  protocol.  Chapter  V 
contains  numerical  results  which  allow  comparisons  among  the 
performances  of  the  protocols  be  made.  Finally,  in  Chapter 
VI,  we  conclude  this  work  and  discussed  some  prospects  for 
future  study.  This  thesis  also  contains  appendices  which 
include  listings  of  key  programs  used  in  computer  simulation 
and  some  numerical  data. 
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II.  STATIC  M-ARY  SPLITTING 


In  this  chapter,  the  study  of  the  static  m-ary  splitting 
for  resolving  collisions  in  a  multichannel  system  occurs.  We 
focus  on  theoretical  analysis  of  this  protocol  in  this 
chapter  and  defer  numerical  examples  together  with 
discussions  until  Chapter  V. 

A.  GENERAL  OPERATION 

In  this  section,  the  m-ary  splitting  algorithm  is 
stated.  We  consider  an  infinite  number  of  independent, 
identical,  bursty  packet  transmitting  users  which 
collectively  form  a  Poisson  traffic  model.  Also  errorless, 
time-slotted  multichannel  is  assumed. 

Let  "t"  be  the  current  time  and  let  stations,  each  of 
which  generated  a  packet  during  the  period  (tQ,tQ+A)  send 
their  packets  at  time  "t".  Each  of  these  stations  randomly 
selects  one  of  the  a  channels  and  sends  its  packet.  Assume 
each  user  knows  whether  his  transmission  is  a  success  or  not 
immediately  by  listening  to  the  channel  after  sending  a 
packet.  When  collision  is  detected,  colliding  users 
involved  in  the  current  slot  of  all  channels  are  divided 
into  "m"  groups  by  flipping  an  m-faced  dice.  The  packets  in 
the  first  group  are  retransmitted  in  the  very  next  slot.  The 
second  group  will  be  transmitted  after  collision,  if  any, 
among  those  in  the  first  group  has  been  cleared,  then  the 
third  group,  and  so  on. 


B.  DERIVATION  OF  Ljj 

Let  Y^j  denote  the  probability  that  j  users  succeed 
when  i  users  attempt  their  retransmissions  in  a  slot  which 
contains  a  cells  (channels).  Of  course  0  <  j  5  i.  At  this 
point,  the  following  question  arises: 

Question:  distribute  i  packets  to  a  channels,  what  is 
the  probability  that  j  of  the  channels  each  contains  exactly 
one  packet  while  each  of  the  others  contain  either  zero  or 
more  than  one  packet  ?  If  a  channel  slot  contains  only  one 
packet,  it  will  be  called  a  success. 

The  answer  can  be  found  as  follows:  let  X°j  denote  the 
number  of  arrangements  in  distributing  i  users  (packets) 
into  a  channels  which  results  in  j  success,  then 
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j  out  of  i  packets,  distributing  j  packets  into  j 
channels  so  that  each  of  them  is  a  success.  The  factor 
(a"j)^"'^  then  tells  the  number  of  arrangements  in 
distributing  the  remaining  i-j  packets  into  the 
remaining  a"j  channels  which  results  in  no  success  at  all. 
Notice  that  is  expressed  in  recursive  form.  Numerical 
values  of  Y,-  z  will  be  discussed  in  Chapter  V. 


of  Binary  Solittin 


Suppose  that  n  packets  have  collided  and  denotes 
the  average  number  of  slots  required  to  solve  this 

collision.  In  this  case,  the  problem  is  equivalent  to  the 
following;  after  a  collision  is  detected,  the  users 
involved  randomly  and  independently  flip  a  fair  coin  to 
determine  whether  to  retransmit  in  the  next  slot  or  wait 
until  later.  Of  course,  in  retransmiting  a  packet,  one  of 
the  a  channels  will  be  picked  randomly. 

•  Remark:  when  a  =  1  *  and  if  no  one  decides  to 

retransmit  immediately  (in  other  words,  if  an  empty 
slot  is  observed),  Massey  s  idea  [Ref.  7]  is  to  restart 
the  process  immediately  in  order  to  avoid  a  guaranteed 
collision.  However,  when  a  >  1  .  it  seems  unnecessary 
to  follow  Massey's  idea. 


Thus 


minca.  i) 


^n  ~  ^  j  Jq  ^  ^n-ij^n-i- j'*’  ^ij^i-j^  (2.2) 

The  2  in  (2.2)  results  from  the  fact  that  each  splitting 
requires  at  least  two  slots.  In  other  words,  for  binary 
splitting,  we  need  at  least  two  slots  to  resolve  a 
collision.  Obviously,  Lq  =  =  0  .  Furthermore,  P(n,i)  is 

defined  to  be 


P(n,i)  =  (1/2)^ 


(2.3) 


which  represents  the  probability  that  i  of  the  n  users 
decide  to  try  again  immediately.  Substitute  (2.3)  into 
(2.2)  we  have 


_  n  rmin(p:,n-i)  minva,'; 

^n-ij^n-i-j"^  ^ij^i-j 


min(o,i) 


(2.4) 


Through  algebraic  manipulation,  we  can  obtain  in 

recursive  form  as  follows: 
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2 .  of  Ternary  Splitting  (m=3) _ 

First  define 

P(n,i;L>^2)  ~  [n! /{ij^!i2!  (n-i]^-i2) !  >]  (1/3)^ 


(2.6) 


Physically  P(n,i^,i2)  represents  the  probability  of  dividing 
n  users  into  three  groups  of  sizes  i^,  i2  and  n-i]^-i2. 
Proceed  similarly  to  binary  splitting,  we  have 


n  n-«, 
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3.  L„  of  General  M-ary  Splittin 


In 

general,  we  define 

P(n; 
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L  =  m  +  2  P(n;i)[  2  2  j  L.  ] 
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From  (2.9),  we  obtain 
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(2.10) 


and  Lq  =  =  0.  Notice  that  (2.10)  is  expressed  in 

recursive  form. 

The  result  of  for  general  m-ary  splitting  can  be 
further  simplified  as  follows: 

—  n-1  „  ,  .  ,  _  .niin(a>i) 
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Finally,  from  (2.11)  we  obtaiji 


min(o,n) 


‘■n  °  i5, 
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(2.12) 


The  form  shown  in  (2.12)  is  very  convenient  for  numerical 
calculations . 

The  following  example  helps  verify  the  correctness 
of  (2.12).  In  this  example  we  consider  a=2,  m=2  and  n=2. 
Figure  2.1  shows  the  state  transition  diagram  for  this 
example.  States  such  as  (1,1)  represents  the  splitting  of 


Figure  2.1  State  Transition  Diagram 
of  Binary  Splitting  with  N=2 


Differentiate  L2(z)  obtained  above  and  set  2  =  1,  we  obtain 
L2  =  [(0.75)(1.5)  >  (0.75)(0.5)3/(0.75)2 
=  8/3 

While  from  (2.12),  we  have 

^2  °  C2  *  2(0.5)2(Y2iLi  •  YjjL,,)  '  0] 
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C.  Ljj 

We  use  L^  to  denote  the  actual  average  number  of  slots 
required  in  successfully  transmitting  a  total  of  n  packets 
including  the  initial  slots  which  contain  collisions  (if 
any ) .  Then 


^  *  j?o 


(2.13) 


This  expression  is  valid  for  all  the  protocols 
considered  in  this  study. 
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D.  LINEAR  BOUNDING  OF  L 


N 


By  linear  bounding,  we  mean  we  assume  there  exist 
^o  ’  i^ln  »  i^un  ^  such  that  satisfies 


U(no-n)Lj^  +  U(n-no  +  l)  ^  ^n 

S  U(ng-n)Lj^  +  U(n-no  +  l)(^_  ^-C) 
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In  practice  it  is  reasonable  to  assume  nQ>a  •  Then 
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Since  U(n„-n+j )+U(n- j -no+1)  =  1,  we  have 
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In  (2.16),  the  terms  which  include  C  can  be  arranged  to 
become 


min  (n, O’) 


..-.  _  j  _  .  min(ot,i) 

md/m)^  2  Y-.+m  2  (.^)  (l/ra)i(l-l/m)"-^  2  Y.  . 

j  =  1  i  =  o^'  j=o 

=  m  -  m(l/m)^  (2.17) 

Similarly,  terms  containing  can  be  combined  to  obtain 


uno 


min(n,tf) 

m(l/m)"  .2  Yj,j(n-j) 


J  =  l 


n-l  _  .  .min(i«a) 

*  m  2  (^(l/m)^(l-l/m)^-^  2  Y.^i-j) 

j=o  ■' 


i  =  0 


(2.18) 


Since 


,2  'Tnj  =  1  -  V„o 


we  have 


^.2^Ynj(n-j)  =  n(l-Yj^g)  -2^  jY^^j 


Let  us  now  obtain  the  closed-form  representation  of 


2  JY 


j=i 


nj 


•  -  ' 

V"  N\% 
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v  v‘.  v  v  j 


The  result  can  be  obtained  by  answering  the  follow  question. 
Given  n  packets  and  a  channels  such  that  nS<x  •  Distribute 
these  packets  into  the  channels  in  a  purely  random  manner. 
What  is  the  expected  number  of  channels,  each  of  them 

a 

containing  exactly  one  packet,  i.e.,  2  ?  Define 

j=o 

=  I  1  if  channel  i  contains  exactly  one  packet 
I  0  otherwise 

then 

P(Xi  =  1)  =  n(l/a)(l-l/a)"-^ 
therefore 

i  jY_.  =  2  E(Xi)  =  n(l-l/a)''‘^  (2.19) 

j=0  1=1 


As  a  matter  of  fact,  this  result  is  good  not  only  for 
n  S  a  ,  but  also  for  n  <  a  ,  since 

2  JY-.  =  an(l/a)(l-l/a)'"’^ 
j=o 

=  nCl-l/a)""'^ 

Thus  (2.18)  now  becomes 

n-1  _  .'"''’(a.i) 

2  (  j )  (l/m)^(l-l/m)^'^  2.  Y^j  (i- j  )  -  ( l/m)"Yj^0^] 

i  =  0  j  =0 

=  n-m(l)"Y„0'^-m  2  (  ■ )  (s)  ^  ( l-li?)''- H  ( 1- 1/a )  ^' ^  (2.20) 

I  =o 

The  second  term  of  (2.20)  can  be  reduced  to 

m  2  (”)  (l/m)^(l-l/m)”'^i(l-l/a)^“^  =  n(l-l/ma)"’^ 

jSQ 

therefore  (2.20)  now  can  be  written  as 


Subsitute  (2.21)  into  (2.15)  and  apply  n  S  n  , 
tedious  algebraic  manipulation,  we  obtain 


^n  ^  t  ™  “  (m-l)C 


m!n(a.n} 


■  .2.  Y^jU(no-n+j  )  (n-j  ) 


‘o 

n-i 


J 


*  m  _2  (;')(l/m)i(l-l/m)^-^  2  YnjU(no-i-j  )  (i-j  ) } 


isO  ' 


0 


minict.n) 

+  {m(l/m)^  2^  Y^jU(no-n+j )  (Lj^_j+  C) 

n-i  -  ,  .  ,  .  minca.  i) 

*  .2(i)(^)^(l-fn)'''^  .2^  YijU(n^-i>j)(Li. -C)}] 


so 


Thus 


A 


un. 


sup  [  A  /  B  ] 

n  in* 


where  A  =  m  -  (m-l)C 


n  „  ...  I  ^  ,  mmw.i) 


1-u  J  “ 0 


0  = 


n  „  ,  .  ,  .  mincflt.i; 

m  2  (i)(-)Ml-f}f)'"-^  2  YijU(no-i-j)(i-j) 

a»o  j  =  0 


Also,  we  have 


^ln„  = 


(2.24) 


inf  [  A  /  B  ] 


Now  we  need  to  select  a  proper  value  for  C.  The 
selection  of  C  can  be  done  as  follows.  For  sufficiently 
large  n,  each  m-ary  splitting  results  in  m  groups.  Because 
collided  users  in  each  group  are  resolved  collectively  and 
each  group  has  approximately  n/m  users,  we  have 

3  m  +  m 


•  claim  :  »  jSn  -  m/(m*l) 

•  check  :  m  +  m  =  m  ♦  m[j3(n/m)-m/ (m-1)] 

=  jSn  -  m/(m-l) 

thus  for  m-ary  splitting,  C  is 
C  =  m/ (m- 1 ) 


(2.25) 


If  we  substitute  (2.25)  into  (2.23),  then  the  first  term  in 
the  numerator  of  (2.23)  will  be  m  -  (m  -  1)C  =  0. 

E.  THROUGHPUT  ANALYSIS 

In  this  section,  maximal  achievable  throughput  X^ax 
obtained.  For  this  purpose,  use  the  results  obtained  in  the 
preceding  section  we  have 


^In  "  ■  C  -  Ln  ^  -  C  V  n>n„>a 


We  can  generalize  (2.26)  to 


Oo-I 

n-C  +  2  C.  5i  V  n>0 
°  i=0  u  n 


(2.26) 


(2.27) 


where  C 


il  =  l-i  -  •  C 


Equation.  (2,13)  can  be  modified  for  this  protocol  as  follows 


rd 


min(a,n) 


0  ^ 

Subsituting  (2.28)  into  (2.27)  gives 
mim.«.n)  no-l 


1  *  .2^  ’TnjtPln  ("-Jl-C  *.2  Ci  Lj 

mincn.a)  n^-i 

j -0  0  ^ 


It  can  be  proven  that 
min(c«.n) 

where  n(l-l/a)”'^  is  only  good  for  a  ^  2.  Also 
rain  (Of, n) 


j 

and 


min(oi,n;  n<,-i 

,?o  1^=0  Vi. -J 


(2.28) 


(2.29) 


(2.30) 


Let  the  mean  arrival  rate  of  the  Poisson  process  be  \  , 
and  N  denote  the  random  variable  representing  the  number  of 
packets  arrived  during  (tQ  ,  tg  +  A)»  we  clearly  have 

E[N]  =  XA 

P[N  =  n]  =  e’^^XA)''  /  n! 


Now  take  the  expected  value  for  all  parts  in  (2.29).  From 
(2.30)  we  have 


V  •  w  ^  -rw  ir»  ir».  ^  nr^.  »r^  k  ^  i  •.  »f*.v^  u  ■  •  ”  i.  ^»ii 


2n(l-l/a)""V^^(\A)''/n! 

n=o 


=  AA  e-^A 

next  we  have 


oo  min(or.n)  n,,-! 

2  [  2  Y^i  2  Ci  5i  „.j]  e-'^“aA)"/n 

n=o  j=o  i»o 


J=w  !»( 

inin(o<.n)  no-i 


ru-i  inini«.nj  no-i 
=  2  [  2  Y^j  2  Ci  5i  n-jl  ® 

nso  j=o  "->  i=o  -J 


■'^^aA)"/n! 
■^^(XA)''/n! 


Ho+O*"'  «  Ma-I  -vA 

>2  [  2  Y„j  2  Ci  5i  n-j]  (XA)"/n! 

n=no  j=n-r\,t«"-'  i^o  J 


Let  us  now  define 


min(o(.n) 

no-l 

LI 

=  J?o  ''"J 

2 

i=o 

Ci^Si 

a 

no-> 

L2 

=  2  Y„j 

2 

i*o 

Ci^<5i 

tnin(jr  n) 

no-t 

U1 

=  Jo 

2 

1=0 

Ci  8i 

(X 

Ho-I 

U2 

•  ^  ^nj 

j  =  n-nofl 

2 

i  =  0 

Ci  5i 

■^U  ^ 

Then  (2.29)  implies 


no- 


1  -  C  >j3in  x(l-e’^/°<)  +2  LI  e'^x^/nl 
o  n=o 

Po-fd-i 

>  2  L2  e-V/nl  <  E[lJ] 

n  -  n« 


^  1  -  C 


n,-i  no+d-> 

♦  2  U1  e'^x"/n!  +  2  U2  e’V/n! 

n=o  n=no 
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Clearly,  it  is  sufficient  that  E[L^]5A  guarantees 

the  stability  of  the  protocol.  On  the  other  hand,  if 
— d 

E[Ljj]>A  )  the  protocol  is  definitely  unstable.  Then  due  to 
the  upper  bound  in  (2.33),  we  also  conclude  that  the 
following  condition  is  sufficient  for  stability: 

no"  I  riotd-i 

+  2  U1  e'*x^/n!+  2  U2  e‘*x^/n!  <  A  (2.34) 

n=o  n^no 

From  (2.33)  and  (2.34),  we  obtain  the  following  in  a 
straightforward  manner.  For  stability  of  the  m-ary 
splitting,  it  is  sufficient  that  the  input  rate  \  satisfies 
the  following  condition: 

no'i 

>  2  U1  e‘^x"/n!  +  2  U2  e'*x^/n!}]  (2.35) 

n  =  0  flario 

Then 

^max  =  sup  [  X  /  ) 

no-1  no+«-i 

+  2  ui  e'V/nl  +  2  U2  e*  V/n!  }]  (2.36) 

n=o  n-no 

The  maximization  of  (2.36)  is  straightforward,  and 
provides  us  with  the  maximum  achievable  throughput,  as  well 
as  the  optimal  A  choice  for  the  m-ary  splitting.  Indeed, 
maxima  is  obtained  at  x  -»  oo  ,  or  equivalently  at  A  -♦  oo 


III.  DYNAMIC  M-ARY  SPLITTING 


In  this  chapter  dynamic  m-ary  splitting  is  considered 
for  resolving  collisions  in  a  multichannel  system.  For 
dynamic  m-ary  splitting,  we  will  assume  that  the  existence 
of  an  infinite  energy  detector.  Again  we  focus  on 
theoretical  analysis  and  defer  numerical  results  until 
Chapter  V. 


DERIVATION  OF  L 


N 


As  mentioned  before,  dynamic  m-ary  splitting  can  be 
analyzed  in  the  same  state  as  static  m-ary  splitting  except 
the  application  of  an  infinite  energy  detector.  Upon  the 
availability  of  an  infinite  energy  detector,  we  can  have 
for  this  protocol  as  follows: 


Lj,  =  Cl/{l-m(n)  (l/m(n)  )^Yj^0^K  n>(n) 


min(Q(.n) 

ra(n)  (l/m(n) )”  2  Y^^jL^^.j 


J  =  l 


n-1 


minW,!; 


II  -  «  .  I  ,  ^ 

*  2^  ’'ijH-J  ]  (3-1) 


where  m  is  a  function  of  n.  One  possible  selection  of  m(n) 
is 


m 


(n)  =  [  n  /  a  ] 


At  this  expression,  [x]  means  the  smallest  integer  >  x. 


B.  LINEAR  BOUNDING  OF  Ljj 

Regarding  the  linear  bounding  of  of  dynamic  m-ary 

splitting  in  (3.1),  they  can  be  obtained  from  (2.23)  by 
subsituting 


m(n)  ->  m 

c(n)  =  m(n)  /  (m(n)  - 1}  -♦  C 
Thus ,  we  have 


sup  [ (m(n) - (m(n) - 1)0 
n  ino 


m(n)  2^  (  I )  (  ^^^^  )^(1-  ^  ^2^  Y^jU(nQ-i  +  j  )  (L^_j+C)  } 


.  miil(d.l) 


/  {n(l-- 


vn-i 


m(n)o( 


.miow  ?) 


+m(n)2^  ( i )  (  m(n)^  ^  ^2^  Y^jU(nQ-i-^  j  )  (i- j  ) }]  (3.2) 


C .  THROUGHPUT  AN ALY  SIS 

In  this  section,  maximal  achievable  throughput 
dynamic  m-ary  splitting  is  obtained.  Using  the  same 
procedure  as  in  Chapter  II,  we  have 

i^lno"  ■  -  ^n  -  Puuo^  ■  n>no>a  (3.3) 

As  far  as  is  concerned  for  dynamic  m-ary  splitting, 

we  can  use  the  same  expression  as  in  (2.28).  Then  we  have 

E[L^]  =  2  P{N  =  n)  (3.4) 


Substitute  (2.28)  into  (3.4)  we  obtain 


E[i4]  S  1  \A(l-e'-^^/®)P{N=n} 


rio-l  _ 

^  2  [  2  Y^j(Ln-j-i3un  (n-j )  )]P{N=n} 

n=o  3=0  o 

no«i'i  «  _ 

n-Mo  j=n-n^+l  o 

no+’<<-l  n-rjo 

-  2  C  2  YniC(n-j)]P{N=n} 

n=no  J=o 

-  2  [2  Y^jC(n-j)]P{N=n) 

n=noM-i  J’O 

Defining,  the  possible  values  of  C(n),  if  n>a  : 
C(n)  >  1 

and  substituting  the  same  expression  of  P[N 
previous  chapters,  we  have 


(3.5) 


=  n]  as  in 


E[L^]  ^  Pun 


rio-l  niin(n,<x) 


•  2 

no+'<-l  oo  _  lA  „ 

*  2  C  2  Y^j(l>L„.j-i3_  (n-j))]  e'^^(XA)^/n! 

j^n-n,+i  jo 


(3.6) 


Similar  to  Chapter  II,  “  A  guarantees  system 

stability.  Thus, 

Ho"!  minifl-n)  _  .. 

nofflt-t  «  _  -iA, 

^2  [2  Y^j(UL^. (AA)"/n!<A  (3.7) 

With  algebraic  manipulation  and  define  x  =  AA  •  Finally,  we 


•■*p  ^ 


=  sup  X  /  CjSun 


rjo-i  fflin(«.n)  _ 

w,?o  ^  .?«  ^nj  ^  *x”/n! 


r»=0  j=o 

no+«<-l  a 


n»no  ,j-n“rt*+'i  "  ‘ 


As  a  matter  of  fact,  by  comparing  (2.36)  and  (3.8) 
observe  that  (3.8)  is  the  result  of  (2.36)  with  C  =  1  . 


IV.  BINARY  SPLITTING  WITH  ENERGY  DETECTOR 


Through  numerical  examples  to  be  given  in  Chapter  V,  we 
shall  see  that  the  best  static  m-ary  splitting  is  binary 
splitting.  Therefore,  in  this  chapter  we  study  the  binary 
splitting  with  energy  detector  protocol.  Only  theoretical 
analysis  of  this  protocol  is  focused  in  this  chapter  and 
numerical  values  and  discussions  are  deferred  until  the  next 
chapter. 

A.  GENERAL  OPERATION 


As  stated  in  Chapter  I,  we  assume  the  availability  of 
infinite  energy  detectors.  After  a  collision  is  detected, 
"n”  users  involved  in  a  collision  divided  into  two  groups  by 
flipping  a  biased  coin.  We  shall  find  the  optimal 
probability  that  a  user  will  join  the  first  group,  which 
minimizes  L^.  To  transmit  a  packet,  one  of  the  available 
channels  will  be  picked  randomly.  The  users  involved  in  the 
second  group,  decide  not  to  retransmit  immediately,  will  be 
held  until  collisions,  if  any,  among  those  in  the  first 
group  are  cleared  and  then  transmit. 


B.  DERIVATION  OF  L 


N 


We  will  use  the  same  Y^j  of  (2.1).  With  this  value  of 


Yj^j  ,  can  be  obtained 


_  n  .  .  min(®.i) 

l-n  =  2  .  ,2  ^2^  YijH.j 

mined,  n-0 

^n-ij^n-i-j^ 

j  =0 


(4.1) 


where  ns2  and  Lq  =  Lj^  =  0.  The  2  in  (4.1)  results  from  the 
same  fact  stated  on  p.  26  .  Equation  (4.1)  can  be  rewritten 
as  follows 


_  ^  mined,  n) 

J=  I 


.  .  mln(a.O 

^2^nCiaA(l-a^)"■"-nC^-ieT^-"(l-(r^)">  .2^  (^-2) 

I  s  1  L  S  O 


J  =  i  ■■  . .  ■•  ■■  j  =  o 

In  (A. 2),  we  assume  that  have  been  optimal. 

If  we  adopt  Massey's  idea  of  MCCRA  [Ref.  8]  in  this 


protocol,  i.e.,  in  each  binary. splitting,  if  the  first  slot 
is  observed  to  be  empty  or  all  users  collide  again,  then  we 
have 


=  [l/{l-YnOa”-(l-an)"}]Cl 
min(d.n) 


n-i  .  . 

^  n*"i‘^n^^'‘^n)^  ^*r\^n~±^n  ^ij^i-j^  (^>3) 

vj  “O 


i3( 


This  method  shall  be  referred  to  as  the  modified  binary 
splitting  with  energy  detector.  As  in  MCCRA,  can  be 
reduced  remarkably  by  using  (4.3).  Regarding  of  L^i  we  can 
use  same  the  expression  as  in  (2.13)  for  this  protocol. 


C.  LINEAR  BOUNDING  OF 

will  consider  the  linear  bounding  of  L^^,  i.e., 

Sin."  -  C  ^  =  Sun."  -  C 

Since  a  >  the  number  of  available  channels,  cannot  be  too 
big  for  practical  interest  we  assume  np>a  ,  and  define 

P(n,i,an)  =  n^i 


1.  No  S  N  S  No  +  o 


For  this  case  we  obtain  the  expression  of  as 

follows : 

j=n-n,i-i 

or-i  «  flo-l  a 

>  2  w  2  YijLi.j  >  2  w  .2  YijLi.j 

i=i  j=o  ^  ^  j*o  ^ 

n-i  i-rto  «  I 

^2  w  [  2  2  YijLi.j]  (4.5) 

i*no  j  =  o  j=;-ndH  ' 


where  W  is  defined  as  follows 


W  =  P(n,i,crn)  P(n,n-i  .oTj^) 


and  substitute  (2.26)  into  (4.5),  we  have 


-  [af .  ( "i" ”y„j  (n- J  )  -  C] 


a-1  ' 


n«-i  a 


>  2  w  2  YijLi.j  >  2  w  2  YijLi.j 

1=1  j  =  o  J  1  =  1  j  =  o  J 

♦  "f'w[’2°Yij(5„n  (i-j)-C).  I  YijLi.j]'  (4.6) 

i:n,  j.o  ‘J  "o  j-i-n,n  ‘  J  ) 

Substitute  the  same  expression  for  ^  Y_:(n-j)  and  2  Y  . 

j=(  “J  j*i 

as  in  Chapter  II  and  apply  (2.19),  after  tedious  algebraic 
manipulation,  we  obtain 


(4.7) 


where  A  =  2  Yj^j(L„_j^C) 

n.-(  mincot.i;  i=n-n,+i 

>  2  W  2  Yij(Li.jH.c) 

i=l  j=0  ^  J 

n-i  "^flC 

+  2  W  2  Y^AL..*C) 

B  =  n[£r„(l-an/a)'''^-(l-a^,){l-(l-(7n)/a}^‘^] 

*  [cT^d-a^)^]  .  2  Y^j(n-J) 
n.-l  iiiin(of.i) 

>  2  W  2  Yidi-j) 

1  =  1  -' 

From  this  expression,  we  obtain  the  upper  bound  and  lower 
bound  of  Lj^, 


un. 


=  sup  [A/B] 

fV>ini  rio+d 

=  inf  [A/B] 

nainirio+i 

N  >  No  ■*•  g 


(4.8) 


Using  defined  W,  for  this  case  can  be  described 
as  follows: 


^n  "  [^/ ( 1- Yj^g ( l-Oj^)  "Yj^qo^}]  .2 


*  [an"dl-an)"l  .2  Y„jL„.j 


a- 1  i 


n,-<  d 


-  .2  W  2  YijLi.j>  2  W  2  YijLi.j 

i  =  |  jso  ^  j  i=(X  j  =  o  J 

I -Mo  « 

•  .2  we  .2  YuLi-j*  2  VijLi.j] 
I -Mo  j-0  J  =  l-no+l 

n-1  «  / 

V  2  w  2  YijLi.j 

i=n,-«(+i  j=o  -'1 


(4.9) 


By  substituting  (2.26)  into  (4.9),  we  obtain 


L„S[l/{l-Y„o(l-cT„)"-Y„oa?)]  }  2 


j 


a-i  « 


Oo-i  oc 


2 

i=t 

rio-t-ot 

i-rto 

2 

W[  2  Yij{;3„n.( 

i=no 

J=0 

n-( 

OC 

2 

W  2  Yi,[3„„  . 

sMo+tf+l  j=0 

♦.  2  W  2  Yij[3„„^(i-J)-C]  j  (4.10) 

i=Ho-rt+l  j=o  ' 

After  some  algebraic  manipulation,  we  obtain  the  upper  bound 


of  Ljj  for  this  case, 


jSun  =  [N  /  D]  ( 

n>no+Q(tl 

no-l  min(a.i) 

where  N  =  1  +  On  (l-CTj,)”  +  2  W  2  Y-.  (Li..:  +  C  ) 

•jsil  j=0  -J 

no+fl[  d 
i=no  jsi-n,+i 


(4.11) 


D  =  nCo^d-an/a)''’^  *  (l-an){l-(l-(rn)/a>''  1 


Ho"  I  min(o(.U 
+  2  W  2  Yi:(i-j) 

i  =  l  j  =  0  -* 


Obviously  lower  bound  of  should  be 


j3in  =  inf  [N  /  D] 


(4.12) 


Now  we  need  to  find  the  value  of  C  for  this 
protocol.  For  a  sufficiently  large  "n" ,  the  "n"  colliding 
users  will  be  divided  approximately  evenly  into  la  groups 
each  with  approximately  n/2a  users,  i.e.. 


S  2  -  2a  L 


n/20( 


V'n  % 


I  « 

■'T^ — =~r« 

••'.v'.'-d 
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claim 


jSn  -  2/(2a-l) 


•  check  :  2  +  2aLj^y2fl(  =  2  +  2a  []8n/2a  -  2/(2a-l)] 

=  j3n  -  2/(2a-l) 


Therefore,  a  better  selection  of  C  for  this  protocol  is 


C  =  2/(2a-l). 


V.  NUMERICAL  RESULTS  AND  DISCUSSION 


We  tackled  the  analysis  of  the  proposed  multichannel 
collision  resolution  algorithm  in  three  protocols  as 
described  in  the  previous  chapter.  We  will  present  the 
corresponding  numerical  results,  compare  with  computer 
simulation  and  discuss  them  accordingly  in  this  chapter. 

There  are  quite  a  number  of  high  level  programming 
languages  available  today.  Some  programming  languages  are 
specifically  developed  for  simulations  of  complex 
organizations  and  systems,  and  others  are  designed  for 
computer  assisted  instruction. 

FORTRAN  is  one  of  the  most  widely  used  languages  for 
expressing  mathematical  relationships.  It  has  excellent 
mathematical  capability  and  is,  therefore,  fully  equipped 
and  supported  by  our  computer  center.  In  this  thesis, 
simulation  is  required  to  verify  the  results  of  theoretical 
analysis.  Therefore,  FORTRAN  is  quite  suitable  for  this 
purpose . 


A.  RESULTS  FOR  M-ARY  SPLITTING 

Through  theoretical  analysis  of  m-ary  splitting,  we  have 
several  expressions  for  L^^,  »  etc.  To  compute  L^  for 
high  enough  value  of  n,  some  numerical  analysis  techniques 
were  needed.  We  obtained  corresponding  to  the  number  of 
available  channel,  as  called  a  •  In  this  thesis,  a  was 
selected  the  value  of  2,  3,  5  and  10.  Numerical  values  of 
Yj^j  ,  corresponding  to  the  above  value  of  a  .  are  provided 
in  Appendix  B.  Some  of  them  are  small  enough  to  be 
neglected,  so  those  are  represented  as  0,  but  the  actual 
value  zero  of  Y^j  is  deleted  from  this  table. 


The  simulation  program,  written  in  FORTRAN,  is  provided 
in  Appendix  A.  For  dividing  groups  and  selecting  channels, 
pseudorandom  number  generator  routine  is  used  with  different 
seed  number  at  each  time.  is  simulated  in  this  program. 

Numerical  values  of  are  compared  with  the  simulation 
results.  The  corresponding  comparisons  are  given  in  Figure 
5.1,  5.2,  5.3  and  5.4  with  each  different  number  of 
available  channels.  As  shown  in  these  figures,  correctness 
of  theoretical  analysis  is  certified,  because  two  results 
are  quite  close  to  each  other.  Three  values  of  ”m”  are 
considered,  i.e.,  2,  3  and  5.  For  each  case  of  a  >  binary 
splitting  shows  the  best  performance  and  ternary  splitting 
is  next  and  so  on.  Also,  it  can  be  seen  from  these  graphs 
that  the  performance  improves  as  the  number  of  available 
channels,  a  is  increased. 

Table  I  contains  the  analytical  results  for  the  linear 
bounding  of  Lj^.  As  shown  in  this  table,  if  we  select  a  much 
smaller  value  of  no  compare  to  n  (number  of  colliding 
users),  the  lower  bound  of  Lj^,  approaches  to  zero.  To 
obtain  the  reasonable  linear  bounding  of  the  differences 
between  those  two  values,  n^  and  n,  should  be  decreased. 

The  analytical  results  for  throughput  are  contained  in 
Table  II.  In  those  papers  mentioned  in  Chapter  I, 
throughput  should  be  less  than  1.0,  because  only  single 
slotted  channel  was  used  in  those  algorithms.  Since  the 
multichannel  random  accessing  problem  is  considered  in  this 
study,  the  results  obtained  here  cannot  be  compared  with  the 
results  in  those  papers.  As  shown  in  Table  II,  throughput 
is  gradually  increased  as  the  number  of  channels  are 
increased.  The  best  performance  was  obtained  with 
parameters  of  a  =  10  and  m  =  2.  Even  with  the  same 
parameters,  different  throughputs  were  obtained  in 
accordance  with  the  value  of  n^. 
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STATIC  M-ARY  SPLITTING 

®  n  :  :  :  ^  : 


n 
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o 
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Figure  5.3  for  M-ary  Splitting  (a  =  5) 
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TABLE  I 

LINEAR  BOUNDING  OF  Lj,  FOR  M-ARY  SPLITTING 


No 

N 

Ajnp 

/^Ino 

10 

150 

1.6027 

0.0 

21 

110 

1.59967 

1.575186 

21 

100 

1.59967 

1.59656 

51 

150 

1.5996387 

1'.  5996025 

10 

150 

1.7144 

0.00015 

15 

150 

1.709172 

0.815373 

21 

110 

1.709029 

1.706522 

21 

150 

1.709029 

1.706522 

10 

150 

2.121764 

1.485665 

15 

150 

2.121758 

2.0847387 

21 

150 

2.1164 

2.0808 

10 

150 

1.2141 

0.00 

21 

150 

1.185685 

0.00 

21 

110 

1.185685 

0.00 

31 

150 

1.18529 

0.00874 

21 

100 

1.185685 

0.018 

21 

60 

1.185685 

1.18056 

51 

150 

1.18528 

1.18523 

10 

150 

1.3564 

0.00 

21 

150 

1.3198 

0.3150 

21 

110 

1.3198 

1.2928 

51 

150 

1.3190 

1.3145 

10 

150 

1.7077 

0.028 

15 

150 

1.7069 

1.3924 

21 

150 

1.6954 

1.6486 

10 

150 

0.9172 

0.0 

10 

30 

0.9172 

0.4112 

21 

50 

0.8250 

0.7794 

41 

150 

0.815959 

0.0134 

51 

150 

0.815897 

0.809036 
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LINEAR  BOUNDING 

TABLE  I 

OF  Lj^  FOR  M-ARY  SPLITTING (CONT'D. 

5 

3 

10 

150 

1.0572 

o 

o 

5 

3 

21 

150 

0.95286 

0.0003 

5 

3 

31 

150 

0.95056 

0.75214 

5 

3 

41 

150 

0.95056 

0.94339 

5 

3 

51 

150 

0.949854 

0.942393 

5 

5 

10 

150 

1.341835 

0.001398 

5 

5 

21 

150 

1.29949 

1.1474 

5 

5 

31 

150 

1.276177 

1.209858 

5 

5 

41 

150 

1.2758 

1.209859 

5 

5 

51 

150 

1.27583 

1.21490 

10 

2 

10 

150 

0.70137 

0.0 

10 

2 

21 

150 

0.54 

0.0 

10 

2 

31 

150 

0.5071 

0.0 

10 

2 

41 

150 

0.4934 

0.0 

10 

2 

51 

150 

0.4886 

0.022 

10 

2 

61 

150 

0.488 

0.4578 

10 

3 

10 

150 

0.8026 

o 

o 

10 

3 

.  21 

150 

0.6684 

0.0 

10 

3 

31 

150 

0.6104 

0.032 

10 

3 

41 

150 

0.6035 

0.5267 

10 

3 

51 

150 

0.6035 

0.5890 

10 

5 

11 

150 

1.0212 

0.00105 

10 

5 

21 

150 

0.9399 

0.4404 

10 

5 

31 

150 

0.8975 

0.7782 

10 

5 

41 

150 

0.8587 

0.7782 

10 

5 

51 

150 

0.85836 

0.77818 
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TABLE  II 

RESULTS  OF  THROUGHPUT  FOR  M-ARY  SPLITTING 


0 

0 

0 

0 

0 


M 

N 

AjOo 

^max 

2 

10 

1.6027 

0.75965567 

2 

21 

1.59967 

0.75965568 

2 

51 

1.5996387 

0.75965568 

3 

10 

1.7144 

0.65249 

3 

15 

1.709972 

0.65249 

5 

10 

2.121764 

0.50916 

2 

10 

1.2141 

0.99682 

2 

21 

1.185685 

0.996823 

3 

10 

1.3564 

0.828139 

3 

21 

1.3198 

0.828130 

5 

10 

1.7077 

0.627832 

5 

15 

1.7069 

0.6278323 

2 

10 

0.9172 

1.39498 

3 

10 

1.0572 

1.1289654 

5 

10 

1.3418 

0.843187 

2 

10 

0.70137 

2.3319745 

3 

10 

0.8026 

1.8328123 

5 

10 

1.012 

1.3495743 

5 

21 

0.9399 

1.3097 

5 

41 

0.8587 

1.3096 

B.  RESULTS  OF  DYNAMIC  M-ARY  SPLITTING 


The  advantage  of  dynamic  m-ary  splitting  is  to  select 
the  "m"  in  m-ary  splitting  as  a  function  of  "n" .  In  other 
words,  static  m-ary  splitting  uses  fixed  m,  while  dynamic 
m-ary  splitting  uses  m(n)  to  minimize  L^. 

Figure  5.5  contains  the  results  of  L^  which  are  obtained 
from  theoretical  analysis  for  this  protocol.  Through  this 
graph,  it  can  also  be  seen  that  the  performance  is  increased 
by  increasing  the  number  of  channel.  As  seen  before,  in  the 
static  m-ary  splitting,  binary  splitting  case  shows  the  best 
performance.  To  compare  with  binary  splitting  and  dynamic 
m-ary  splitting,  both  data  were  plotted  in  same  graph  as 
shown  in  Figure  5.6.  In  Figure  5.6,  dynamic  m-ary  splitting 
outperforms  static  binary  splitting. 

Results  of  linear  bounding  of  L^^  for  this  protocol  are 
given  in  Table  III.  Also  in  throughput  analysis,  dynamic 
ra-ary  splitting  performs  better  than  static  m-ary  splitting. 
These  results  are  shown  in  Table  IV. 


Figure  5.5  Results  of  for  Dynamic  M-ary  Splitting 
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TABLE  III 


LINEAR 

BOUNDING 

OF  Ljj 

FOR  DYNAMIC  M-ARY 

SPLITTING 

Pf 

la 

N 

^  1  no 

2 

10 

150 

1.272472 

1.178883 

2 

21 

150 

1.2117635 

1.178896 

2 

31 

150 

1.202447 

1.178896 

3 

10 

150 

0.951344 

0.831984 

3 

21 

150 

0.873781 

0.832144 

3 

31 

150 

0.858532 

0.832144 

5 

■  10 

150 

0.745224 

0.559788 

5 

21 

150 

0.615479 

0.562131 

5 

31 

150 

0.592985 

0.56213176 

10 

11 

150 

0.54172959 

0.3127287 

10 

21 

150 

0.410154 

0.330307 

10 

31 

150 

0.375021 

0.330331 

10 


41 


150 


0.359558 


0.330331 


TABLE  IV 

RESULTS  OF  THROUGHPUT  FOR' DYNAMIC  M-ARY  SPLITTING 


a 

N 

^max 

1 

10 

1.272472 

0.855374 

1 

21 

1.217635 

0.855386 

3 

10 

0.951344 

1.2111 

3 

21 

0.873781 

1.2114 

5 

10 

0.745224 

1.7848 

5 

21 

0.615479 

1.7998 

5 

31 

0.592985 

1.7998 

10 

11 

0.54172959 

2.9477998 

10 

21 

0.410154 

3.0952285 

10 

31 

0.375021 

3.097495 

10 


41 


0.359558 


3.0974956 


C.  RESULTS  OF  BINARY  SPLITTING  WITH  ENERGY  DETECTOR 


In  accordance  with  the  availability  of  an  infinite 
energy  detector,  if  collision  occurs,  the  total  number  of 
users  involved  in  the  collision  are  known  exactly.  With 
this  advantage,  the  optimal  which  minimizes  Lj^,  were 
sought  for  each  case  of  a  and  n.  For  most  of  the 
cases,  is  very  close  to  0.5.  If  we  use  expression  in 
(4.2),  which  is  expressed  in  recursive  form,  with  (7^  =  0.5, 
the  results  of  L^^  are  exactly  the  same  as  in  static  binary 
splitting.  With  this  result  in  mind,  simulation  is  not 
needed  for  this  protocol.  The  results  of  L^  for  binary 
splitting  with  energy  detector  are  given  in  Figure  5.7. 

If  we  use  modified  binary  splitting  with  energy  detector 
protocol,  as  expressed  in  (4.3),  value  of  L^^  can  be  reduced 
remarkably.  Figure  5.8  shows  these  results.  However,  this 
modified  protocol  is  very  sensitive  to  channel  errors,  as 
in  MCCRA,  because  the  decision  of  retransmission  depends  on 
the  feedback  information  strongly. 

Analytic  results  of  L^  are  plotted  in  Figure  5.9,  5.10, 
5.11  and  5.12  for  several  values  of  a  •  L^,  the  actual 
average  number  of  slots  required  in  successfully 
transmitting  a  total  of  n  packets  including  the  initial 
slot,  is  commonly  used  for  all  protocols  considered  in  this 
thesis . 


BINARY  SPLITTING 


Figure  5.7  for  Binary  Splitting  with  Energy  Detector 
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MODIFIED  BINARY  SPLITTING 


AVERAGE  NUMBER  OF  SLOTS 


VI.  CONCLUSION 


This  thesis  has  presented  a  few  collision  resolution 
algorithms  for  a  multichannel  random  accessing 
communications  system.  Since  some  multi channeling  can  be  a 
feasible  configuration,  it's  worth  the  effort  to  observe  the 
application  of  CRA  to  a  multichannel  system.  It  has  been 
noted  that  there  are  a  few  ways  to  handle  the  multichannel 
CRA.  In  this  thesis  we  focus  on  one  class  of  methods  which 
resolve  collisions  in  all  channels  jointly. 

We  consider  three  different  collision  resolution 
algorithms.  These  algorithms  have  been  analyzed  to  obtain 
the  desired  system  performance  by  computing  L^,  linear 
bounding  of  and  maximal  achievable  throughput.  Through 
the  theoretical  analysis  and  computer  simulations,  we  have 
observed  that  binary  splitting  shows  the  best  performance  in 
static  m-ary  splitting.  However,  in  applying  infinite 
energy  detectors  to  m-ary  splitting,  we  are  able  to  select 
the  proper  value  of  m,  as  a  function  of  a  and  n,  which 
minimizes  the  and  maximizes  the  throughput  for  each  case. 
Thus,  we  have  algorithm  better  than  the  best  static  m-ary 
splitting.  This  better  algorithm  is  dynamic  m-ary  splitting. 

Since  it  was  proven  that  the  best  static  m-ary  splitting 
is  binary  splitting,  the  binary  splitting  with  energy 
detector  was  also  considered  and  evaluated.  The  unique 
feature  of  this  algorithm  is  the  optimal  probability  (t^ 
that  a  user  will  join  the  first  group  when  the  collision 
involves  "n”  users  occurs.  Modified  binary  splitting  with 
energy  detector  was  also  analyzed  by  adopting  Massey's  idea. 
We  observed  that  this  modified  algorithm  outperforms  the 
dynamic  m-ary  splitting  except  in  the  system  stability.  If 
we  consider  only  the  ideal  channel,  certainly  this  modified 


algorithm  will  be  the  best  one.  In  practice,  however,  it  is 
impossible  to  expect  an  ideal  channel. 

When  considering  channel  noise  or  channel  errors,  (as 
far  as  we  are  concerned  in  the  multichannel  collision 
resolution  algorithms),  we  conclude  that  dynamic  m-ary 
splitting  will  maintain  minimum  and  maximum  achievable 
throughput.  We  note  again  that  the  throughputs  computed  here 
cannot  be  compared  with  those  in  several  previous  papers 
[Refs.  6,  7,  8,  9,  10,  12],  because  the  multichannel  system 
is  much  different  from  the  single  channel  system.  These 
algorithms  were  shown  to  be  stable  -  to  have  maximum 
throughput  within  a  certain  arrival  rate. 

As  mentioned  earlier,  the  network  performance  can  be 
evaluated  by  three  properties  -  throughput,  delay  and 
stability.  The  maximum  achievable  throughput  and  the  system 
stability  were  found,  but  not  for  average  delay,  in  this 
thesis.  Therefore,  these  proposed  algorithms  can  be  further 
studied  for  average  delay  in  the  future.  However,  the 
results  developed  here  are  applicable  to  a  general  system  in 
which  a  central  facility  is  accessible  by  a  number  of 
independent  users.  If,  for  example,  the  number  of  users  that 
the  facility  (multichannel  system)  can  handle  simultaneously 
is  less  than  the  maximum  number  of  users  that  can  place 
requests  upon  it,  then  collisions  will  arise  that  may  be 
resolved  with  the  techniques  developed  here. 


APPENDIX  A 

SIMULATION  PROGRAM  FOR  EVALUATING  M-ARY  SPLITTING  ALGORITHM 


DEFINITION  OF  VARIABLES  ************************ 


A  =  NUMBER  OF  AVAIABLE  CHANNELS 
N  =  NUMBER  OF  COLLIDING  USERS 

M  =  NUMBER  OF  SPLITTING 

LN  =  TOTAL  LENGTH  OF  SLOTS  REQUIRED  TO 
,  RESOLVE  COLLISION 
L(N)=  AVERAGE  LENGTH  OF  SLOTS 
.T  =  NUMBER  OF  COMPUTING  TIMES 


.DECLARATION  OF  VARIABLES 


********************** 


INTEGER  N.I,T,A.T1,T2 
INTEGER  T:3.T4,T5.t6,M 
REAL  L(0: 100) ,C,LN,TD,SC 

SELECT  THE  NUMBER  OF  CHANNEL  AND  SPLITTING 


A  =  2 ,  3 ,  5  or  10 
M  =  2’  3,  or  5 
T  =  2000  OR  MORE 
L(0)  =  0. 

l1i)  =  1.  ^ 

TD  =  FLOAT  if  T) 

DO  31  N  =  2,30 
LN  =  0. 

DO  41  I  =  1,T 
C  =  0. 

GO  TO  33 
C  =  C  +  1. 

33  IF(M.EQ.2)  THEN 

CALL  TIL\CH2(A,N,I,C,T1,T2) 

GO  TO  34 

ELSE  IFfM.EQ. 3)  THEN 

CALL  TRACH3(a,N,I,C,T1,T2,T3) 

GO  TO  35 
ELSE 

CALL  TRACH5(A,N,I,C,T1,T2,T3,T4,T5) 
GO  TO  36 
,  END  IF 

34  IFffTl.EQ.NKOR.  (T2.EQ.N)  )  GO  TO  30 

if(a.eq.2)  then 

CALL  C0NSA2(A,N,I,C,T1,L,SC) 

c  =  c  +  sc 

CALL  CONSA2(A,N,I,C,T2,L,SC) 

c=c+sc 

ELSE  IF(A.EQ.3)  THEN 

CALL  CONSA3(a,N,I,C,T1,L,SC) 

c  =  c  +  sc 

call  C0NSA3(A,N,I,C,T2,L,SC) 

c  s  c  +  sc 

ELSE  IF(A.EQ.5)  THEN 

CALL  CONSA5(a,N,I,C,T1,L,SC) 

c  =  c  sc 

CALL  CONSA5(A,N,I,C,T2,L,SC) 

C  =  C  +  SC 
ELSE 


CALL  C0NSA0(A,N,I,C,T1,L,SC) 

c  —  c  sc 

CALL  CONSAO(A,N,I,C,T2,L,SC) 
C  =  C  +  SC 
END  IF 
GO  TO  42 


^  CALL  XONSA2?A,N,I,C,Tl,L,SC) 

C  =  C  +  SC 

CALL  CONSA2(A,N,I,C,T2,L,SC) 

C  =  C  +  SC 

CALL  CONSA2(A,N,I,C,T3,L,SC) 

C  =  C  +  SC 

ELSE  IF(A.EQ.3)  THEN 

CALL  C0NSA3(A,N,I,C,T1,L,SC) 

C  =  C  +  SC 

CALL  C0NSA3(A,N,I,C,T2,L,SC) 

C  =  C  +  SC 

CALL  CONSA3(A,N,I,C,T3,L»SC) 

c  —  c  sc 

ELSE  IF(A.EQ.5)  THEN 

CALL  CONSA5(A,N,I,C,Tl,L,SC) 

c  =  c  +  sc 

CALL  CONSA5(A,N,I,C,T2,L,SC) 

C  =  C  +  SC 

CALL  CONSA5(A,N,I,C,T3,L,SC) 

C  =  C  +  SC 
ELSE 

CALL  CONSAO(A,N,I,C,T1,L,SC) 

c  —  c  sc 

CALL  CONSAO(A,N,I,C,T2,L,SC) 

c  —  c  +  sc 

CALL  CONSAO(A,N,I,C,T3,L,SC) 

C  =  C  +  SC 
END  IF 
GO. TO  42 

IF((T1.EQ.N) .OR. (T2.EQ.N) .OR. (T3.EQ, 
iT4.Eg.N) .OR. (T5.EQ.N))  GO  TO  30 
IF(A.EQ.2j  THEN 

CALL  C0NSA2(A,N,I,C,T1,L,SC) 

C  =  C  +  SC 

CALL  C0NSA2(A,N,I,C,T2,L,SC) 

C  =  C  +  SC 

CALL  C0NSA2(A,N,I,C,T3,L,SC) 

C  =  C  +  SC 

CALL  C0NSA2(A,N,I,C,T4,L,SC) 

C  =  C  +  SC 

CALL  C0NSA2(A,N,I,C,T5,L,SC) 

c  —  c  sc 

ELSE  IFfA.EQ.S)  THEN 

CALL  C0NSA3(A,N,I,C,T1,L,SC) 

c  =  c  +  sc 

CALL  C0NSA3(A,N,I,C,T2,L,SC) 

c  =  c  +  sc 

CALL  C0NSA3(A,N,I,C,T3,L,SC) 

c  —  c  sc 

CALL  C0NSA3(A,N,I,C,T4,L,SC) 

c  s  c  +  sc 

CALL  CONSA3(A,N,I,C,T5,L,SC) 

c  =  c  +  sc  '  ’  ’  ’  ’  ’  ’  ’ 

ELSE  IF(A.EQ.5)  THEN 

CALL  C0NSA5(A,N,I,C,T1,L,SC) 

c  =  c  +  sc 

CALL  C0NSA5(A,N,I,C,T2,L,SC) 

C  =  C  +  SC 

CALL  CONSA5(A,N,I,C,T3,L,SC) 

C  =  C  +  SC 


(T3.EQ.N))  GO  TO  30 
SC) 


N) .OR. 


CALL  C0NSA5(A,N,I,C,T4,L,SC) 

c  “  c  sc 

CALL  CONSA5(A,N,I,C,T5,L,SC) 

C  =  C  +  SC 
ELSE 

CALL  CONSAO(A,N,I,C,T1,L,SC) 

C  =  C  +  SC 

CALL  CONSAO(A,N,I,C,T2,L,SC) 

c  —  c  +  sc 

CALL  CONSAO(A,N,I,C,T3,L,SC) 

C  s  Q  + 

CALL  CONSAO(A,N,I,C,T4,L,SC) 

c  —  c  sc 

CALL  CONSAO(A,N,I,C,T5,L,SC) 

C  =  C  +  SC 
END  IF 
GO  TO  42 

WHEN  ALL  USERS  COLLIDED  AGAIN,  RESTART  SPLITTING 


C  =  C  +  FLOAT(M-l) 

TC  =  N 

IF(A.EQ.2)  THEN 

CALL  C0NSA2 ( A . N , I , 
IF(SC.EQ.O)  GO  to 


C^TC,L,SC) 


IF(SC.EQ.O)  GO  TO  32 

c  =  c  +  sc 

ELSE  IF(A.EQ.3)  THEN 

CALL  C0NSA3(A.N.I,C.TC,L,SC) 

ifTsc.eq.O)  Gd  to  3i 
c  —  c  +  sc 

ELSE  IF(A.EQ.5)  THEN 

CALL  CONSA5(A.N,I,C.TC,L,SC) 

IF(SC.EQ.0rG6  to  3i 

C  =  C  +  SC 

ELSE 

CALL  CONSAO(A,N.I,C.TC,L,SC) 

IF(SC.EQ.O)  Gd  to 
C  =  C  +  SC 

END  IF 

2  LN  =  LN  +  C 

1  CONTINUE 

COMPUTE  Ljj  AND  PRINT  IT  OUT 

L(N)  =  LN  /  TD 
WRltE(4,100)N,L(N) 

1  CONTINUE 

00  F0RMAT(4X,I3,3X,F10.4) 

STOP 
END 

SUBROUTINE  OF  BINARY  SPLITTING  **************«>->* 


SUBROUTINE  TRACH2 ( A , N , I , C , TCI , TC2 ) 
INTEGER  J.TC1,TC2;N,I  A 
REAL  NRUdOKd 
REAL*8  DSEED 

DSEED  =  182669. pO+DFLOAT(I)*12 

+DBLE(C)"14152+DFL0AT (A)*328 
CALL  GGUBS (DSEED, N, NR) 

TCI  =  0 
TC2  =  0 

DO  201  J  =  l.N 

IFtNR(J) .tT.0.5)  THEN 
tci  =  TCI  +  1 


ELSE 
TC2  = 
END  IF 
CONTINUE 
RETURN 


TC2  +  1 


END 


SUBROUTINE  OF  TERNARY  SPLITTING  ****•***■****•***■*■*•* 


SUBROUTINE  TRACH3 ( A , N , I , C , TCI , TC2 , TC3 ) 
INTEGER  I. J,TC1,TC2,TC3,N,A 
REAL  NRfldoKc 


REAL*8  DSEED 
DSEED  =  2240 
+DBLE 


7.D0+DFLOAT(I)*108 
__^_^C]“13502tDFLOAT(A)*490 
CALL  GGUBSXdSEED,N,NR) 

TCI  =  0 
TC2  =  0 
TC3  =  0 

DO  301  J  =  l.N 

IFiNR(J) .LT. 0.3334)  THEN 
TCI  =  TCI  +  1 

ELSE^IF(NRj[j).LT.  0.6667)  THEN 


301 


TC2  = 
ELSE 
TC3  = 
END  IF 
CONTINUE 
RETURN 
END 


TC^ 
TC3 


SUBROUTINE  OF  SPLITTING  WITH  FIVE  BRANCHES 


«.'« 


SUBROUTINE  TRACH5 lA,N,I,C,TCl,TC2,TC3,TC4,TC5) 
INTEGER  I . J , TC 1 , TC2 , Tfi3 , TtA , TfiS , N  M , A 
REAL  NR(i6o5,C 

T5t?AT  -nOCTA’ 


if 


REAL*8  DSEED 
DSEED  =  198463. DO+DFLOATfI)*327 

+  DBLE  f  C)*8  9  5  0  +  DFLOAT ( A ) *8  3 1 
CALL  GGUBS( DSEED, N, NR) 

TCl  =  0 
TC2  =  0 
TC3  =  0 
TC4  =  0 
TC5  =  0 

DO  501  J  =  l.N 

IFrNR(J) .LT.0.2)  THEN 
TCI  =  TCI  +  1 
ELSE  IFCNRfj) .LT.0.4) 

TC2  =  TC2  +  1 
ELSE  IF(NRTj) .LT,0.6) 


501 


Tn  =  'Tr 
ELSE  IF(NR?J) 
TC4  =  TC4  + 
ELSE 

TC5  =  TC5  + 
END  IF 
CONTINUE 
RETURN 
END 


,  LT . 0 . 8 ) 
1 


THEN 

THEN 

THEN 


DISTRIBUTE  COLLIDING  USERS  INTO  10  CHANNELS  * 

# 


# 


SUBROUTINE  ASHIO ( A , N , I , C . TCI . TC2 . TC3 , TC4 , TC5 , 
TC6 , TC7 , TC8 , TC9 , TClO ) 

INTEGER  I , TCI , TC2 , TC3 , TC4 , TC5 , TC6 , TC7 , TC8 , J 
INTEGER  T(!:9,t610,A,N 
REAL  NR(100j,C 
REAL*8  DSEED 

DSEED  =  123719.DO+DFLOAT(I)*119 

+  DBLErc)*9171+ DFLOAT ( A ) *5  0  8 
CALL  GGUBS (DSEED, N, NR) 

TCI  =  0 
TC2  =  0 
TC3  =  0  • 


THEN 

THEN 

THEN 

THEN 

THEN 

THEN 

THEN 

THEN 


TC4  =  0 
TC5  =  0 
TC6  =  0 
TC7  =  0 
TC8  =  0 
TC9  =  0 
TCIO  =  0 
DO  601  J  =  l.N 

IFrNR(J) .LE.0.1)  THEN 
TCI  =  TCI  +  1 

ELSE  IF(NRrJ) .LE.0.2)  THEN 
TC2  =  TC2  +1  ^ 

ELSE  IF(NR7j) .LE.0.3)  THEN 
TC3  =  TC3  +1 

ELSE  IF(NRTj)  .LE.0-.4)  THEN 
TC4  =  TC4  t  1 

ELSE  IF(NR(J).LE.0.5)  THEN 
TC5  =  TC5  t  1 

ELSE  IF(NR7j) .LE.0.6)  THEN 
TC6  =  TC6  +  1 

ELSE  IF(NR(J) .LE.0.7)  THEN 
TC7  =  TC7  +  1 

ELSE  IF(NRij) .LE.0.8)  THEN 
TC8  =  TC8  +  1 

ELSE  IF(NRTj) .LE.0.9)  THEN 
TC9  =  TC9  +  1 
ELSE 

TCIO  =  TCIO  +  1 
END  IF 

01  CONTINUE 
RETURN 
END 

DISTRIBUTE  COLLIDING  U?5R? . MP  SECCEEDING  USERS , 
IN  TWO  CHANNELS  **'*'WW**'*'********WW**'s’f  *********** 

SUBROUTINE  CONS A2  (A , N , I , C . T , L , S C ) 

INTEGER  N,T,CH1.CH2  I,A  cdLt 
REAL  C,SC,L( 0:100) 

IF(T.L^:.lj  GO  TO  433 

CALL  TRACH2(A,T,I.C,CH1.CH2) 

IFUCHl.EQ.N)  .OR.  (CH2.E(i.N))  GO  TO  43 
COLL  =  0 

IFTcHI.LE.1)  GO  TO  12 
COLL  =  COLL  +  CHI 

2  IF(CH2.LE.l)  GO  TO  13 

COLL  =  cOlL  +  CH2 
SC  =  1,+  L(COLL) 

GO  TO  44 

3  SC  =  1,  +  L(COLL) 

GO  TO  44 

3  SC  =  0. 

GO  TO  44 

33  SC  =  1. 

4  RETURN 
END 

DISTRIBUTE  COLLIDING  USERS  AND  SUCCEEDING  USERS 
IN  THREE  CHANNELS  «''f*****'********'*Vf**i!f******Vf*** 

SUBROUTINE  C0NSA3 (A , N , I , C , T , L . SC ) 

INTEGER  N , A , T , I . CHI , CH2  Cft3 , cOlL 
REAL  C,SC,l[0: 100) 

IF(T.Lfe.lj  Go  TO  455 
call  TRACH3(A,T,I,C,CH1.CH2,CH3) 
IF(TcHI.EQ.N) .OA. (cft2.E0.N).OR. (CH3.EQ.N)) 
y/  GO  TO  45 

COLL  =  0 

IF(CHl.LE.l)  GO  TO  14 
COLL  =  cOlL  +  CHI 


GO  TO  43 


14 

IF(CH2.LE.l 
COLL  =  O 

15 

IF(CH3.LE.l 
COLL  =  Cl 
SC  =  1,+ 
GO  TO  46 

16 

SC  =  1.+  L(i 
GO  TO  46 

45 

SC  =  0, 

GO  TO  46 

455 

SC  =  1. 

46 

RETURN 

END 

GO  TO  15 
.L  +  CH2 
GO  TO  16 
L  +  CH3 
JCOLL) 


DISTRIBUTE  COLLIDING  USERS  AND  SncCEEDING  USERS 
,IN  FIVE  channels 

SUBROUTINE  CONSA5 ( A , N , I , C , T . L , SC ) 

INTEGER  N , I . T , CHI , CH2 , CH3 , Cft4 , CHS , A , COLL 
REAL  C,SC,L?0:100) 

IF(T.LE.1J  GO  TO  477 

CALL  TRACKS  U , T , I , C , CHI . CH2 , CH3 , CH4 , CHS ) 
IFTTcHI . EQ . N ) . or . ( CH2 . EQ . N ) ! OR . iCH3 . EQ . N ) . OR 
y/  (CH4.EQ.N)  .OR.  (CH5.EQ.N))  GO  TO  47 

IF(CHi.LE.l)  GO  TO  21 
COLL  =  COLL  +  CHI 

11  IF(CH2.LE.l}  GO  TO  22 

COLL  =  COLL  +  CH2 

12  IF(CH3.LE.l}  GO  TO  23 

COLL  =  COLL  +  CH3 
Z3  IF(CH4.LE.l)  GO  TO  24 

COLL  =  COLL  +  CH4 
lU  IF(CH5.LE.l)  GO  TO  25 

COLL  =  COLL  +  CHS 
SC  =  1.+  L(COLL) 

GO  TO  26 

15  SC  =  1.+  L(COLL) 

GO  TO  26 
+77  SC  =  1. 

GO  TO  26 
+7  SC  =  0. 

Z6  RETURN 

,  END 

DISTRIBUTE  SUCCEEDING  U$ERS  AND  COLLIDING  USERS 
IN  TEN  CHANNELS  ^  ^  ^  ^  ^  ^  .ww ^  ^  ^ +v »» 

SUBROUTINE  CONSAO (A , N , I , C , T , L, SC ) 

INTEGER  N , I , T , A , CHI . c62 ! cft3  Cft4 , CHS , CH6 , CH7 
INTEGER  efts ,Cft9, CHI 0,COtL 
REAL  C.SC.LfO: lOO) 

IF(T.LE.ll  GO  TO  18 

CALL  ASHlO ( A , T , I , C , CHI , CH2 . CH3 , CH4 , CHS , CH6 , 

#  _  cft7’cft8!CH9!CHl6) 

IF( (CHl.EQ.N) .0R.TCH2.feQ.N} .OR. (CH3 . EQ. N) . OR 

#  (CH4.EQ.N) .OR. fCHS.EQ.N) .OR. fCH6 .EQ.N) .OR 

#  fCH7.EQ.N; .OR.fCHS.EQ.NJ.OR. (CH9 . EQ. N) . OR 

#  (CHIO.EQ.N))  GO  TO  19 


CCH7.EQ.N 
(CHIO.EQ. 
COLL  =  0 
IF(CHl.LE.l) 
COLL  =  C6 


EQ.N) .OR 
EQ . N ) . OR 
EQ.N) .OR 


IF(CHl.LE.l)  GO  TO  1 
COLL  =  COLL  +  CHI 
IF(CH2.LE.l}  GO  TO  2 
COLL  =  COLL  +  CH2 
IF(CH3.LE.l)  GO  TO  3 
COLL  =  COLL  +  CH3 
IF(CH4.LE.l)  GO  TO  4 
COLL  =  cOlL  +  CH4 
IF(CHS.LE.ll  GO  TO  S 
COLL  =  cOlL  +  CHS 


IF(CH6.LE.l 
COLL  =  C( 
IF(CH7.LE.l, 


COLL  =  COLL 


GO  TO  6 
ILL  +  CH6 
GO  TO  7 


IF(CH8.LE.l)  GO 
COLL  =  COLL  + 
IF(CH9.LE.l}  GO 
COLL  =  COLL  + 
IF(CH10.LE. ij  GO 
COLL  =  COLL  + 


COLL  = 
SC  =  ] 
GO  TO 
SC  =  1.+ 
GO  TO  11 
SC  =  1. 

GO  TO  11 
SC  =  0. 
RETURN 
END 


GO  TO  8 
LL  +  CH8 
GO  TO  9 
LL  +  CH9 
T  GO  TO  10 


ilcOLL) 


L(COLLj 


APPENDIX  B 

NUMERICAL  VALUES  OF  Yjj 


Y^j  denotes  the  probability  that  j  users  succeed  when  i 
users  attempt  their  retransmissions  in  a  slot  which 
contains  a  channels.  Using  the  equation  (2.1),  we  computed 
Yjj^j  corresponding  to  each  value  of  a  (1  to  10)  where  n 
(collising  users)  is  up  to  150.  Since  the  total  number  of 
these  values  is  too  big  to  attach  here,  this  appendix 
contains  only  some  of  them  for  the  value  of  a  is  2,  3  and 
5.  As  has  been  noted  before,  some  of  the  values  are  small 
enough  to  be  neglected,  so  they  are  represented  as  0  in 
following  table.  On  the  other  hand,  some  probabilities  are 
very  close  to  1.0,  Hence  they  are  represented  as  1.0,  These 
values  are  valid  for  all  protocols  studied  in  this  thesis. 


TABLE  V 

Yjj  (a  =  2) 


i 

j_ 

a 

i 

j 

Ill 

T 

1 

O.lOOOOOOOD+01 

T 

24 

F 

0.999997140+00 

2 

0 

0.49999997D+00 

2 

24 

1 

0.286102290-05 

2 

2 

0.50000003D+00 

2 

25 

0 

0.999998510+00 

3 

0 

0.25000000D+00 

2 

25 

1 

0.149011610-05 

3 

1 

0.75000000D+00 

2 

26 

0 

0.999999230+00 

4 

0 

0.50000000D+00 

2 

26 

1 

0.774860380-06 

4 

1 

0.50000000D+00 

2 

27 

0 

0.999999600+00 

5 

0 

0.68750000D+00 

2 

27 

1 

0.402331350-06 

5 

1 

0.31250000D+00 

2 

28 

0 

0.999999790+00 

6 

0 

0.81250000D+00 

2 

28 

1 

0.208616260-06 

6 

1 

0. 18750000D+00 

2 

29 

0 

0.999999890+00 

7 

0 

0.89062500D+00 

2 

29 

1 

0.108033420-06 

7 

1 

0. 10937500D+00 

2 

30 

0 

0.999999940+00 

8 

0 

0.937500000+00 

2 

30 

1 

0.558793540-07 

8 

1 

0.62500000D-01 

2 

31 

0 

0.999999970+00 

9 

0 

0.96484375D+00 

2 

31 

1 

0.288710000-07 

9 

1 

0.35156250D-01 

2 

32 

0 

0.999999990+00 

10 

0 

0.98046875D+00 

2 

32 

1 

0.149011610-07 

10 

1 

0. 19531250D-01 

2 

33 

0 

0.999999990+00 

11 

0 

0.98925781D+00 

2 

33 

1 

0.768341120-08 

11 

1 

0. 10742187D-01 

2 

34 

0 

0.100000000+01 

12 

0 

0.99414062D+00 

2 

34 

1 

0.395812090-08 

12 

1 

0.585937500-02 

2 

35 

0 

0.100000000+01 

13 

0 

0.996826170+00 

2 

35 

1 

0.203726810-08 

13 

1 

0.317382810-02 

2 

36 

0 

0.100000000+01 

14 

0 

0.998291020+00 

2 

36 

1 

0.104773790-08 

14 

1 

0.170898440-02 

2 

37 

0 

0.100000000+01 

15 

0 

0.999084470+00 

2 

37 

1 

0.538420860-09 

15 

1 

0.915527340-03 

2 

38 

0 

0.100000000+01 

16 

0 

0.999511720+00 

2 

38 

1 

0.276486390-09 

16 

1 

0.488281250-03 

2 

39 

0 

0.100000000+01 

17 

0 

0.999740600+00 

2 

39 

1 

0.141881170-09 

17 

1 

0.259399410-03 

2 

40 

0 

0.100000000+01 

18 

0 

0.999862670+00 

2 

40 

1 

0.727595760-10 

18 

1 

0.137329100-03 

2 

41 

0 

0.100000000+01 

19 

0 

0.999927520+00 

2 

41 

1 

0.372892830-10 

19 

1 

0.724792480-04 

2 

42 

0 

0.100000000+01 

20 

0 

0.999961850+00 

2 

42 

1 

0.190993890-10 

20 

1 

0.381469730-04 

2 

43 

0 

0.100000000+01 

21 

0 

0.999979970+00 

2 

43 

1 

0.977706800-11 

21 

1 

0.200271610-04 

2 

44 

0 

0.100000000+01 

22 

0 

0.999989510+00 

2 

44 

1 

0.500222090-11 

22 

1 

0.104904170-04 

2 

45 

0 

0.100000000+01 

23 

0 

0.999994520+00 

2 

45 

1 

0.255795380-11 

23 

1 

0.548362730-05 

2 

46 

0 

0. 100000000+01 

TABLE 

V 

Yij  (  a  =  2) 

(CONT'D 

.) 

2 

46 

1 

0.13073986D-11 

2 

71 

0 

0.100000000+01 

< 

2 

47 

0 

0. lOOOOOOOD+01 

2 

71 

1 

0.601393390-19 

V 

w  -  ^ 

2 

47 

1 

0.66791017D-12 

2 

72 

0 

0.100000000+01 

2 

48 

0 

0. lOOOOOOOD+01 

2 

72 

1 

0.304931860-19 

2 

48 

1 

0.34106051D-12 

2 

73 

0 

0.100000000+01 

2 

49 

0 

0. lOOOOOOOD+01 

2 

73 

1 

0.154583510-19 

yAis 

2 

49 

1 

0. 17408297D-12 

2 

74 

0 

0.100000000+01 

4 

2 

50 

0 

0. lOOOOOOOD+01 

2 

74 

1 

0.783505480-20 

v/vj; 

2 

50 

1 

0.88817842D-13 

2 

75 

0 

0.100000000+01 

2 

51 

0 

0. lOOOOOOOD+01 

2 

75 

1 

0.397046690-20 

2 

51 

1 

0.45297099D-13 

2 

76 

0 

0.100000000+01 

2 

52 

0 

0. lOOOOOOOD+01 

2 

76 

1 

0.201170320-20 

2 

52 

1 

0.23092639D-13 

2 

77 

0 

0.100000000+01 

2 

53 

0 

0. ioboooooD+01 

2 

77 

1 

0.101908650-20 

2 

53 

1 

0. 11768364D-13 

2 

78 

0 

0.100000000+01 

2 

54 

0 

0. lOOOOOOOD+01 

2 

78 

1 

0.516160700-21 

;  . 

2 

54 

1 

0.59952043D-14 

2 

79 

0 

0.100000000+01 

•  ■-  - 

2 

55 

0 

0. lOOOOOOOD+01 

2 

79 

1 

0.261389070-21 

2 

55 

1 

0.30531133D-14 

2 

80 

0 

0.100000000+01 

.  J 

2 

56 

0 

0. lOOOOOOOD+01 

2 

80 

1 

0.132348900-21 

2 

56 

1 

0. 15543122D-14 

2 

81 

0 

0.100000000+01 

2 

57 

0 

0. lOOOOOOOD+01 

2 

81 

1 

0.670016300-22 

2 

57 

1 

0.79103391D-15 

2 

82 

0 

0.100000000+01 

2 

58 

0 

0. lOOOOOOOD+01 

2 

82 

1 

0.339144050-22 

2 

58 

1 

0.40245585D-15 

2 

83 

0 

0.100000000+01 

2 

59 

0 

0-.  lOOOOOOOD+01 

2 

83 

1 

0.171639980-22 

2 

59 

1 

0.204697370-15 

2 

84 

0 

0.100000000+01 

2 

60 

0 

0. 100000000+01 

2 

84 

1 

0.868539640-23 

2 

60 

1 

0.104083410-15 

2 

85 

0 

0.100000000+01 

•**  ■* . 

2 

61 

0 

0.100000000+01 

2 

85 

1 

0.439439700-23 

2 

61 

1 

0.529090660-16 

2 

86 

0 

0.100000000+01 

2 

62 

0 

0.100000000+01 

2 

86 

1 

0.222304790-23 

2 

62 

1 

0.268882140-16 

2 

87 

0 

0. 100000000+01 

2 

63 

0 

0 . 100000000+01 

2 

87 

1 

0.112444860-23 

.  '♦  .*-■ . 

2 

63 

1 

0.136609470-16 

2 

88 

0 

0.100000000+01 

2 

64 

0 

0. 100000000+01 

2 

88 

1 

0.568686670-24 

2 

64 

1 

0.693889390-17 

2 

89 

0 

0.100000000+01 

2 

65 

0 

0. 100000000+01 

2 

89 

1 

0.287574510-24 

2 

65 

1 

0.352365710-17 

2 

90 

0 

0.100000000+01 

2 

66 

0 

0.100000000+01 

2 

90 

1 

0.145402840-24 

2 

66 

1 

0.178893360-17 

2 

91 

0 

0.100000000+01 

2 

67 

0 

0. 100000000+01 

2 

91 

1 

0.735092150-25 

2 

67 

1 

0.908019320-18 

2 

92 

0 

0.100000000+01 

2 

68 

0 

0 . 100000000+01 

2 

92 

1 

0 . 371585040-25 

2 

68 

1 

0.460785920-18 

2 

93 

0 

0.100000000+01 

'■.*  *.*•  * 

2 

69 

0 

0. 100000000+01 

2 

93 

1 

0.187812000-25 

2 

69 

1 

0.233781090-18 

2 

94 

0 

0.100000000+01 

‘  *•"  *'«“ 

2 

70 

0 

0. 100000000+01 

2 

94 

1 

0.949157440-26 

**•  / 

2 

70 

1 

0.118584610-18 

2 

95 

0 

0.100000000+01 

*  *  *  ^  • 

TABLE  V 


Yij  (a  =  2) 

(CONT'D 

.) 

2 

95 

1 

0.47962743D-26 

2 

120 

0 

2 

96 

0 

0. lOOOOOOOD+01 

2 

120 

1 

2 

96 

1 

0.24233807D-26 

2 

121 

0 

2 

97 

0 

0.  lOOOOOOOD+01 

2 

121 

1 

2 

97 

1 

0. 1224312 lD-26 

2 

122 

0 

2 

98 

0 

O.lOOOOOOOD+01 

2 

122 

1 

2 

98 

1 

0.61846695D-27 

2 

123 

0 

2 

99 

0 

0. lOOOOOOOD+01 

2 

123 

1 

2 

99 

1 

0.31238892D-27 

2 

124 

0 

2 

100 

0 

O.lOOOOOOOD+01 

2 

124 

1 

2 

100 

1 

0.15777218D-27 

2 

125 

0 

2 

101 

0 

0.  lOOOOOOOD+01 

2 

125 

1 

2 

101 

1 

0. 79674951D-28 

2 

126 

0 

2 

102 

0 

0. lOOOOOOOD+01 

2 

126 

1 

2 

102 

1 

0.40231906D-28 

2 

127 

0 

2 

103 

0 

0. lOOOOOOOD+01 

2 

127 

1 

2 

103 

1 

0.203131680-28 

2 

128 

0 

2 

104 

0 

0. 100000000+01 

2 

128 

1 

2 

104 

1 

0.102551920-28 

2 

129 

0 

2 

105 

0 

0.100000000+01 

2 

129 

1 

2 

105 

1 

0.517689970-29 

2 

130 

0 

2 

106 

0 

0.100000000+01 

2 

130 

1 

2 

106 

1 

0.261310170-29 

2 

131 

0 

2 

107 

0 

0.100000000+01 

2 

131 

1 

2 

107 

1 

0.131887680-29 

2 

132 

0 

2 

108 

0 

0.100000000+01 

2 

132 

1 

2 

108 

1 

0.665601390-30 

2 

133 

0 

2 

109 

0 

0.100000000+01 

2 

133 

1 

2 

109 

1 

0.335882180-30 

2 

134 

0 

2 

110 

0 

0.100000000+01 

2 

134 

1 

2 

110 

1 

0.169481840-30 

2 

135 

0 

2 

111 

0 

0.100000000+01 

2 

135 

1 

2 

111 

1 

0.855112900-31 

2 

136 

0 

2 

112 

0 

0.100000000+01 

2 

136 

1 

2 

112 

1 

0.431408310-31 

2 

137 

0 

2 

113 

0 

0.100000000+01 

2 

137 

1 

2 

113 

1 

0.217630080-31 

2 

138 

0 

2 

114 

0 

0.100000000+01 

2 

138 

1 

2 

114 

1 

0.109778010-31 

2 

139 

0 

2 

115 

0 

0.100000000+01 

2 

139 

1 

2 

115 

1 

0.553704860-32 

2 

140 

0 

2 

116 

0 

0.100000000+01 

2 

140 

1 

2 

116 

1 

0.279259840-32 

2 

141 

0 

2 

117 

0 

0.100000000+01 

2 

141 

1 

2 

117 

1 

0.140833630-32 

2 

142 

0 

2 

118 

0 

0.100000000+01 

2 

142 

1 

2 

118 

1 

0.710186670-33 

2 

143 

0 

2 

119 

0 

0.100000000+01 

2 

143 

1 

2 

119 

1 

0.358102600-33  , 

2 

144 

0 

O.lOOOOOOOD+01 
0. 18055593D-33 
O.lOOOOOOOD+01 
0.91030283D-34 
0. lOOOOOOOD+01 
0.45891299D-34 
0. lOOOOOOOD+01 
0.23133729D-34 
0. lOOOOOOOD+01 
0. 11660904D-34 
0. lOOOOOOOD+01 
0.58774718D-35 
0. lOOOOOOOD+01 
0.29622458D-35 
0. lOOOOOOOD+01 
0. 14928778D-35 
0.  lOOOOOOOD+01 
0.75231638D-36 
0.  lOOOOOOOD+01 
0.37909693D-36 
0.  lOOOOOOOD+01 
0. 19101783D-36 
0. lOOOOOOOD+01 
0.96243600D-37 
0. lOOOOOOOD+01 
0.48489142D-37 
0.  lOOOOO'OOD+01 
0.24428242D-37 
0. lOOOOOOOD+01 
0. 12305956D-37 
0. lOOOOOOOD+01 
0. 61988960D-38 
0.  lOOOOOOOD  +  01 
0.31224069D-38 
0.  lOOOOOOOD  +  01 
0. 15726829D-38 
0.  lOOOOOOOD+01 
0.79208115D-39 
0.  lOOOOOOOD  +  01 
0.39891044D-39 
0.  lOOOOOOOD+01 
0.20089015D-39 
0.  lOOOOOOOD+01 
0. 10116254D-39 
0.  lOOOOOOOD+01 
0.50940002D-40 
0. lOOOOOOOD+01 
0.25649367D-40 
0.  lOOOQOOOD+01 


81 
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Yij  ( a  =  2)  (CONT'D. ) 


2  144  1  G. 12914367D-40 
2  145  0  0. lOOOOOOOD+01 
2  145  1  0.65020249D-41 
2  146  0  0. lOOOOOOOD+01 
2  146  1  0.32734332D-41 
2  147  0  O.lOOOOOOOD+01 
2  147  1  0. 16479270D-41 
2  148  0  0. lOOOOOOOD+01 
2  148  1  0.82956869D-42 
2  149  0  0. lOOOOOOOD+01 
2  149  1  0.41758694D-42 
2  150  0  0. lOOOOOOOD+01 
2  150  1  0.21019477D-42 


TABLE  VI 
Yu  (a  =  3) 


1 

j 

lii 

a 

1 

j 

lii 

T 

1 

G. lOOOOOOOD+Ol 

T 

Ti 

0 

0.974124860+00 

2 

0 

0.33333329D+00 

3 

17 

1 

0.258688230-01 

2 

2 

0.66666671D+00 

3 

17 

2 

-0.631871590-05 

3 

0 

0.11111124D+00 

3 

18 

0 

0.981733100+00 

3 

1 

0 . 66666663D+00 

3 

18 

1 

0.182645270-01 

3 

3 

0.22222213D+00 

3 

18 

2 

0.236951850-05 

4 

0 

0.25925926D+00 

3 

19 

0 

0.987144730+00 

4 

1 

0.29629630D+00 

3 

19 

1 

0.128543850-01 

4 

2 

0.44444444D+00 

3 

19 

2 

0.882761780-06 

5 

0 

0.25925926D+00 

3 

20 

0 

0.990978470+00 

5 

1 

0.493827160+00 

3 

20 

1 

0.902120590-02 

5 

2 

0.24691358D+00 

3 

20 

2 

0.326948810-06 

6 

0 

0.33333333D+00 

3 

21 

0 

0.993684820+00 

6 

1 

0.543209880+00 

3 

21 

1 

0.631506090-02 

6 

2 

0.123456790+00 

3 

21 

2 

0.120454820-06 

7 

0 

0.443072700+00 

3 

22 

0 

0.995589360+00 

7 

1 

0.499314130+00 

3 

22 

1 

0.441059870-02 

7 

2 

0.576131690-01 

3 

22 

2 

0.441667690-07 

8 

0 

0.557384550+00 

3 

23 

0 

0.996925900+00 

8 

1 

0.417009600+00 

3 

23 

1 

0.307408290-02 

a 

2 

0.256058530-01 

3 

23 

2 

0.161243760-07 

9 

0 

0.659807960+00  ' 

3 

24 

0 

0.997861490+00 

9 

1 

0.329218110+00 

3 

24 

1 

0.213850320-02 

9 

2 

0.109739370-01 

3 

24 

2 

0.586340940-08 

10 

0 

0.744449520+00 

3 

25 

0 

0.998514920+00 

10 

1 

0.250978000+00 

3 

25 

1 

0.148507560-02 

10 

2 

0.457247370-02 

3 

25 

2 

0.212442370-08 

11 

0 

0.811106030+00 

3 

26 

0 

0.998970350+00 

11 

1 

0.187031110+00 

3 

26 

1 

0.102965380-02 

11 

2 

0.186285970-02 

3 

26 

2 

0.767153000-09 

12 

0 

0.862012900+00 

3 

27 

0 

0.999287160+00 

12 

1 

0.137241950+00 

3 

27 

1 

0.712837750-03 

12 

2 

0.745143860-03 

3 

27 

2 

0.276175080-09 

13 

0 

0.900098040+00 

3 

28 

0 

0.999507170+00 

13 

1 

0.996084230-01 

3 

28 

1 

0.492826280-03 

13 

2 

0.293541520-03 

3 

28 

2 

0.991397730-10 

14 

0 

0.928178920+00 

3 

29 

0 

0.999659720+00 

14 

1 

0.717069250-01 

3 

29 

1 

0.340284880-03 

14 

2 

0.114155040-03 

3 

29 

2 

0.354944860-10 

15 

0 

0.948661590+00 

3 

30 

0 

0.999765320+00 

15 

1 

0.512944990-01 

3 

30 

1 

0.234679250-03 

15 

2 

0.439057830-04 

3 

30 

2 

0.126766020-10 

16 

0 

0.963478190+00 

3 

31 

0 

0.999838330+00 

16 

1 

0.365050800-01 

3 

31 

1 

0.161667940-03 

16 

2 

0.167260130-04 

3 

31 

2 

0.451695020-11 

3 

32 

0 

TABLE 
Yij  (a  =  3) 
0.99988874D+00 

VI 

(CONT'D 
3  48 

.) 

1 

0.25406922D-06 

3 

32 

1 

0. 11125536D-03 

3 

48 

2 

0.84847710D-19 

3 

32 

2 

0. 16060268D-11 

3 

49 

0 

0.99999983D+00 

3 

33 

0 

0.99992351D+00 

3 

49 

1 

0.17290822D-06 

3 

33 

1 

0.76488060D-04 

3 

49 

2 

0.29486083D-19 

3 

33 

2 

0.56988046D-12 

3 

50 

0 

0.99999988D+00 

3 

34 

0 

0.99994746D+00 

3 

50 

1 

0. 11762464D-06 

3 

34 

1 

0.52537253D-04 

3 

50 

2 

0. 10238223D-19 

3 

34 

2 

0.20183266D-12 

3 

51 

0 

0.99999992D+00 

3 

35 

0 

0.99996395D+00 

3 

51 

1 

0.79984756D-07 

3 

35 

1 

0.36054978D-04 

3 

51 

2 

0.35520367D-20 

3 

35 

2 

0.71354982D-13 

3 

52 

0 

0.99999995D+00 

3 

36 

0 

0.99997528D+00 

3 

52 

1 

0.54368723D-07 

3 

36 

1 

0.24723414D-04 

3 

52 

2 

0. 12313727D-20 

3 

36 

2 

0.25184111D-13 

3 

53 

0 

0.99999996D+00 

3 

37 

0 

0.99998306D+00 

3 

53 

1 

0.36942850D-07 

3 

37 

1 

0. 16940117D-04 

3 

53 

2 

0.42655395D-21 

3 

37 

2 

0.88744011D-14 

3 

54 

0 

0.99999997D+00 

3 

38 

0 

0.99998840D+00 

3 

54 

1 

0.25093257D-07 

3 

38 

1 

0. 11598638D-04 

3 

54 

2 

0. 14765329D-21 

3 

38 

2 

0.31224745D-14 

3 

55 

0 

0.99999998D+00 

3 

39 

0 

0.99999206D+00 

3 

55 

1 

0. 17038631D-07 

3 

39 

1 

0.79359105D-05 

3 

55 

2 

0.51075038D-22 

3 

39 

2 

0. 10970856D-14 

3 

56 

0 

0.99999999D+00 

3 

40 

0 

0. 9999945 7D+00 

3 

56 

1 

0. 11565616D-07 

3 

40 

1 

0.54262636D-05 

3 

56 

2 

0. 17655569D-22 

3 

40 

2 

0.38494232D-15 

3 

57 

0 

0.99999999D+00 

3 

41 

0 

0.99999629D+00 

3 

57 

1 

0.78480968D-08 

3 

41 

1 

0. 37079468D-05 

3 

57 

2 

0.60991964D-23 

3 

41 

2 

0. 13489432D-15 

3 

58 

0 

0.99999999D+00 

3 

42 

0 

0. 99999747D+00 

3 

58 

1 

0. 53238551D-08 

3 

42 

1 

0.25322564D-05 

3 

58 

2 

0.21056749D-23 

3 

42 

2 

0. 47213012D-16 

3 

59 

0 

0. lOOOOOOOD+01 

3 

43 

0 

0.99999827D+00 

3 

59 

1 

0.36104305D-08 

3 

43 

1 

0. 17283654D-05 

3 

59 

2 

0.72651943D-24 

3 

43 

2 

0. 16505362D-16 

3 

60 

0 

0. lOOOOOOOD+01 

3 

44 

0 

0.99999882D+00 

3 

60 

1 

0.24477495D-08 

3 

44 

1 

0. 11790400D-05 

3 

60 

2 

0.25052394D-24 

3 

44 

2 

0.57637771D-17 

3 

61 

0 

0. lOOOOOOOD+01 

3 

45 

0 

0.99999920D+00 

3 

61 

1 

0. 16590302D-08 

3 

45 

1 

0.80389091D-06 

3 

61 

2 

0.86338759D-25 

3 

45 

2 

0.20106199D-17 

3 

62 

0 

0. lOOOOOOOD+01 

3 

46 

0 

0.99999945D+00 

3 

62 

1 

0. 11241516D-08 

3 

46 

1 

0.54783677D-06 

3 

62 

2 

0.29738906D-25 

3 

46 

2 

0.70067058D-18 

3 

63 

0 

0. lOOOOOOOD+01 

3 

47 

0 

0. 99999963D+00 

3 

63 

1 

0. 76152206D-09 

3 

47 

1 

0.37316417D-06 

3 

63 

2 

0. 10237984D-25 

3 

47 

2 

0.24393717D-18 

3 

64 

0 

0. lOOOOOOOD+01 

3 

43 

0 

0.99999975D+00 

3 

64 

1 

0.51573981D-09 

TABLE  VI 
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3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


64 

2 

0.35227472D-26 

3 

81 

0 

0.100000000+01 

65 

0 

0. lOOOOOOOD+01 

3 

81 

1 

0.662497580-12 

65 

1 

0.34919883D-09 

3 

81 

2 

0.438404120-34 

65 

2 

0. 12115268D-26 

3 

82 

0 

0.100000000+01 

66 

0 

0. lOOOOOOOD+01 

3 

82 

1 

0.447117710-12 

66 

1 

0.23638074D-09 

3 

82 

2 

0.149788070-34 

66 

2 

0.41646234D-27 

3 

83 

0 

0.100000000+01 

67 

0 

0. lOOOOOOOD+01 

3 

83 

1 

0.301713570-12 

67 

1 

0. 15997485D-09 

3 

83 

2 

0.511621820-35 

67 

2 

0.14309219D-27 

3 

84 

0 

0.100000000+01 

68 

0 

0. lOOOOOOOD+01 

3 

84 

1 

0.203565790-12 

63 

1 

0. 10824169D-09 

3 

84 

2 

0.174700130-35 

68 

2 

0.491427720-28 

3 

85 

0 

0.100000000+01 

69 

0 

0. lOOOOOOOD+01 

3 

85 

1 

0.137326120-12 

69 

1 

0.73222318D-10 

3 

85 

2 

0.596365920-36 

69 

2 

0. 16869907D-28 

3 

86 

0 

0.100000000+01 

70 

0 

0. lOOOOOOOD+01 

3 

86 

1 

0.926278180-13 

70 

1 

0.49522341D-10 

3 

86 

2 

0.203521700-36 

70 

2 

0.57886935D-29 

3 

87 

0 

0.100000000+01 

71 

0 

0. lOOOOOOOD+01 

3 

87 

1 

0.624699230-13 

71 

1 

0.33486535D-10 

3 

87 

2 

0.694368160-37 

71 

2 

0. 19854939D-29 

3 

88 

0 

0.100000000+01 

72 

0 

0. lOOOOOOOD+01 

3 

88 

1 

0.421253120-13 

72 

1 

0.226387840-10 

3 

88 

2 

0.236838750-37 

72 

2 

0.680740770-30 

3 

89 

0 

0.100000000+01 

73 

0 

0.100000000+01 

3 

89 

1 

0.284026730-13 

73 

1 

0.153021410-10 

3 

89 

2 

0.807611070-38 

73 

2 

0.233305520-30 

3 

90 

0 

0.100000000+01 

74 

0 

0.100000000+01 

3 

90 

1 

0.191478690-13 

74 

1 

0.103411730-10 

3 

90 

2 

0.275321960-38 

74 

2 

0.799287430-31 

3 

91 

0 

0.100000000+01 

75 

0 

0.100000000+01 

3 

91 

1 

0.129070820-13 

75 

1 

0.698727910-11 

3 

91 

2 

0.938363220-39 

75 

2 

0.273728570-31 

3 

92 

0 

0.100000000+01 

76 

0 

0.100000000+01 

3 

92 

1 

0.869927890-14 

76 

1 

0.472029520-11 

3 

92 

2 

0.319738580-39 

76 

2 

0.937088810-32 

3 

93 

0 

0.100000000+01 

77 

0 

0.100000000+01 

3 

93 

1 

0.586255750-14 

77 

1 

0.318826960-11 

3 

93 

2 

0.108921930-39 

77 

2 

0.320692610-32 

3 

94 

0 

0. 100000000+01 

78 

0 

0.100000000+01 

3 

94 

1 

0.395039720-14 

78 

1 

0.215311710-11 

3 

94 

2 

0.370966010-40 

78 

2 

0.109710630-32 

3 

95 

0 

0.100000000+01 

79 

0 

0.100000000+01 

3 

95 

1 

0.266161510-14 

79 

1 

0.145381410-11 

3 

95 

2 

0.126314590-40 

79 

2 

0.375200860-33 

3 

96 

0 

0.100000000+01 

80 

0 

0.100000000+01 

3 

96 

1 

0.179308810-14 

80 

1 

0.981477890-12 

3 

96 

2 

0.430007110-41 

80 

2 

0.128273800-33 

3 

97 

0 

0.100000000+01 
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3 

97 

1 

0. 12078441D-14 

3 

113 

2 

0.46207869D-49 

3 

97 

2 

0.14635330D-41 

3 

114 

0 

0.  lOOOOOOOD+01 

3 

98 

0 

0. lOOOOOOOD+01 

3 

114 

1 

0. 1440763 lD-17 

3 

98 

1 

0.81353070D-15 

3 

114 

2 

0. 15677670D-49 

3 

98 

2 

0.49800775D-42 

3 

115 

0 

0.  lOOOOOOOD+01 

3 

99 

0 

0. lOOOOOOOD+01 

3 

115 

1 

0.96893425D-18 

3 

99 

1 

0.54788802D-15 

3 

115 

2 

0.53183835D-50 

3 

99 

2 

0.16942532D-42 

3 

116 

0 

0. lOOOOOOOD+01 

3 

100 

0 

O.lOOOOOOOD+01 

3 

116 

1 

0.65157318D-18 

3 

100 

1 

0.36894816D-15 

3 

116 

2 

0. 18038962D-50 

3 

100 

2 

0.57627659D-43 

3 

117 

0 

O.lOOOOOOOD+01 

3 

101 

0 

0. lOOOOOOOD+01 

3 

117 

1 

0.43812679D-18 

3 

101 

1 

0.24842510D-15 

3 

117 

2 

0.61175609D-51 

3 

101 

2 

0. 19597285D-43 

3 

118 

0 

0. lOOOOOOOD+01 

3 

102 

0 

0. lOOOOOOOD+01 

3 

118 

1 

0.29458098D-18 

3 

102 

1 

0. 16725650D-15 

3 

118 

2 

0.20743454D-51 

3 

102 

2 

0.66630768D-44 

3 

119 

0 

0. lOOOOOOOD+01 

3 

103 

0 

0. lOOOOOOOD+01 

3 

119 

1 

0. 19805162D-18 

3 

103 

1 

0. 11259751D-15 

3 

119 

2 

0.703268090-52 

3 

103 

2 

0.22650063D-44 

3 

120 

0 

0. 100000000+01 

3 

104 

0 

0. lOOOOOOOD+01 

3 

120 

1 

0.133143940-18 

3 

104 

1 

0.75793796D-16 

3 

120 

2 

0.238395960-52 

3 

104 

2 

0.76980606D-45 

3 

121 

0 

0.100000000+01 

3 

105 

0 

0. lOOOOOOOD+01 

3 

121 

1 

0.895023180-19 

3 

105 

1 

0.51015055D-16 

3 

121 

2 

0.808008720-53 

3 

105 

2 

0.26158458D-45 

3 

122 

0 

0.100000000+01 

3 

106 

0 

0. lOOOOOOOD+01 

3 

122 

1 

0.601613380-19 

3 

106 

1 

0.34333942D-16 

3 

122 

2 

0.273825180-53 

3 

106 

2 

0.88871686D-46 

3 

123 

0 

0.100000000+01 

3 

107 

0 

0. lOOOOOOOD+01 

3 

123 

1 

0.404363090-19 

3 

107 

1 

0.23105231D-16 

3 

123 

2 

0.927837380-54 

3 

107 

2 

0.30188160D-46 

3 

124 

0 

0.100000000+01 

3 

108 

0 

0. lOOOOOOOD+01 

3 

124 

1 

0.271767060-19 

3 

108 

1 

0. 15547445D-16 

3 

124 

2 

0.314349280-54 

3 

108 

2 

0. 10252583D-46 

3 

125 

0 

0.100000000+01 

3 

109 

0 

0. lOOOOOOOD+01 

3 

125 

1 

0.182639160-19 

3 

109 

1 

0. 10460935D-16 

3 

125 

2 

0.106486880-54 

3 

109 

2 

0.34814066D-47 

3 

126 

0 

0.100000000+01 

3 

110 

0 

0. lOOOOOOOD+01 

3 

126 

1 

0.122733510-19 

3 

110 

1 

0. 70379382D-17 

3 

126 

2 

0.360681380-55 

3 

110 

2 

0. 11819590D-47 

3 

127 

0 

0.100000000+01 

3 

111 

0 

0. lOOOOOOOD+01 

3 

127 

1 

0.824717250-20 

3 

111 

1 

0.47346130D-17 

3 

127 

2 

0.122150760-55 

3 

111 

2 

0.40121545D-48 

3 

128 

0 

0.100000000+01 

3 

112 

0 

0. lOOOOOOOD+01 

3 

128 

1 

0.554140730-20 

3 

112 

1 

0.31848448D-17 

3 

128 

2 

0.413632200-56 

3 

112 

2 

0. 13617009D-48 

3 

129 

0 

0.100000000+01 

3 

113 

0 

0. lOOOOOOOD+01 

3 

129 

1 

0.372313300-20 

3 

113 

1 

0.21421872D-17 

3 

129 

2 

0.140048700-56 
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3 

130 

0 

0. lOOOOOOOD+01 

3 

140 

2 

0.93172646D-62 

3 

130 

1 

0.25013297D-20 

3 

141 

0 

O.lOOOOOOOD+01 

3 

130 

2 

0.47412320D-57 

3 

141 

1 

0.31364823D-22 

3 

131 

0 

0. lOOOOOOOD+01 

3 

141 

2 

0.31504420D-62 

3 

131 

1 

0.16803804D-20 

3 

142 

0 

0. lOOOOOOOD+01 

3 

131 

2 

0. 16049131D-57 

-  3 

142 

1 

0.21058179D-22 

3 

132 

0 

O.lOOOOOOOD+01 

3 

142 

2 

0. 10651494D-62 

3 

132 

1 

0. 11288052D-20 

3 

143 

0 

0. lOOOOOOOD+01 

3 

132 

2 

0.54320137D-58 

3 

143 

1 

0.14137651D-22 

3 

133 

0 

0. lOOOOOOOD+01 

3 

143 

2 

0.36008598D-63 

3 

133 

1 

0.75823783D-21 

3 

144 

0 

0. lOOOOOOOD+01 

3 

133 

2 

0. 18383151D-58 

3 

144 

1 

0.94910104D-23 

3 

134 

0 

0. lOOOOOOOD+01 

3 

144 

2 

0. 12171920D-63 

3 

134 

1 

0.50929257D-21 

3 

145 

0 

0. idooooooD+01 

3 

134 

2 

0.62205611D-59 

3 

145 

1 

0.63712801D-23 

3 

135 

0 

0. lOOOOOOOD+01 

3 

145 

2 

0.41140524D-64 

3 

135 

1 

0.34206218D-21 

3 

146 

0 

0, lOOOOOOOD+01 

3 

135 

2 

0.21047011D-59 

3 

146 

1 

0.42768133D-23 

3 

136 

0 

0. lOOOOOOOD+01 

3 

146 

2 

0. 13903973D--64 

3 

136 

1 

0.22973065D-21 

3 

147 

0 

0. lOOOOOOOD+01 

3 

136 

2 

0.71203819D-60 

3 

147 

1 

0.28707377D-23 

3 

137 

0 

0. lOOOOOOOD+01 

3 

147 

2 

0.46985841D-65 

3 

137 

1 

0.15427990D-21 

3 

148 

0 

0. lOOOOOOOD+01 

3 

137 

2 

0.24086230D-60 

3 

148 

1 

0.19268444D-23 

3 

138 

0 

0. lOOOOOOOD+01 

3 

148 

2 

0. 15876494D-65 

3 

138 

1 

0. 10360402D-21 

3 

149 

0 

0. lOOOOOOOD+01 

3 

138 

2 

0.81468132D-61 

■  3 

149 

1 

0. 1293242 4D-23 

3 

139 

0 

0. lOOOOOOOD+01 

3 

149 

2 

0.53641669D-66 

3 

139 

1 

0. 69569847D-22 

3 

150 

0 

0. lOOOOOOOD+01 

3 

139 

2 

0.27552483D-61 

3 

150 

1 

0.86794791D-24 

3 

140 

0 

0. lOOOOOOOD+01 

3 

150 

2 

0. 18122186D-66 

3 

140 

1 

0.46713567D-22 
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i 

j 

Zli 

a 

i 

j 

hi 

T 

0 

O.OOOOOOOOD+OO 

5 

Tl 

0 

0.21120913D+00 

1 

1 

0. lOOOOOOOD+01 

5 

11 

1 

0.44674827D+00 

2 

0 

0.20000000D+00 

5 

11 

2 

0.29257112D+00 

2 

1 

O.OOOOOOOOD+OO 

5 

11 

3 

0.48660477D-01 

2 

2 

0.80000000D+00 

5 

11 

4 

0.81100800D-03 

3 

0 

0.40000219D-01 

5 

12 

0 

0.26103017D+00 

3 

1 

0.48000000D+00 

5 

12 

1 

0.47434335D+00 

3 

2 

O.OOOOOOOOD+OO 

5 

12 

2 

0.23767399D+00 

3 

3 

0.47999978D+00 

5 

12 

3 

0.26709195D-01 

4 

0 

0. 10399974D+00 

5 

12 

4 

0.24330240D-03 

4 

1 

0. 12800045D+00 

5 

13 

0 

0.31870761D+00 

4 

2 

0.57599989D+00 

5 

13 

1 

0.48348588D+00 

4 

3 

O.OOOOOOOOD+OO 

5 

13 

2 

0. 18362252D+00 

4 

4 

0. 19199992D+00 

5 

13 

3 

0.14113701D-01 

5 

0 

0.65600254D-01 

5 

13 

4 

0.70287360D-04 

5 

1 

0.31999977D+00 

5 

14 

0 

0.38150301D+00 

5 

2 

0. 19200021D+00 

5 

14 

1 

0.47462477D+00 

5 

3 

0.38399978D+00 

5 

14 

2 

0. 13660296D+00 

5 

4 

O.OOOOOOOOD+OO 

5 

14 

3 

0.72495657D-02 

5 

5 

0.38399986D-01 

5 

14 

4 

0. 19680461D-04 

6 

0 

0.89919848D-01 

5 

15 

0 

0.44634408D+00 

6 

1 

0.23808024D+00 

5 

15 

1 

0.45125828D+00 

6 

2 

0.40319998D+00 

5 

15 

2 

0.98749596D-01 

6 

3 

0. 15359999D+00 

5 

15 

3 

0.36426742D-02 

6 

4 

0. 11519995D+00 

5 

15 

4 

0.53673984D-05 

7 

0 

0.95039994D-01 

5 

16 

0 

0.51047719D+00 

7 

1 

0.29747204D+00 

5 

16 

1 

0. 41789669D+00 

7 

2 

0.3386&798D+00 

5 

16 

2 

0. 69826522D-01 

7 

3 

0.21503999D+00 

5 

16 

3 

0. 17981609D-02 

7 

4 

0.53760000D-01 

5 

16 

4 

0. 14313062D-05 

3 

0 

0. 11362560D+00 

5 

17 

0 

0.57177287D+00 

8 

1 

0.32727043D+00 

5 

17 

1 

0.37882223D+00 

8 

2 

0.34836478D+00 

5 

17 

2 

0.48529849D-01 

3 

3 

0. 18923519D+00 

5 

17 

3 

0.87467591D-03 

8 

4 

0.21504000D-01 

5 

17 

4 

0.37434163D-06 

9 

0 

0. 13885697D+00 

5 

18 

0 

0. 62878595D+00 

9 

1 

0.36200450D+00 

5 

18 

1 

0.33752450D+00 

9 

2 

0.35721214D+00 

5 

13 

2 

0.33269271D-01 

9 

3 

0. 13418495D+00 

5 

18 

3 

0.42017813D-03 

9 

4 

0.77414400D-02 

5 

18 

4 

0.96259277D-07 

10 

0 

0. 17069313D+00 

5 

19 

0 

0. 68068236D+00 

10 

1 

0.40564738D+00 

5 

19 

1 

0.29656144D+00 

10 

2 

0.33702910D+00 

5 

19 

2 

0.22556512D-01 

10 

3 

0.34049915D-01 

5 

19 

3 

0. 19966996D-03 

10 

4 

0.25804800D-02 

5 

19 

4 

0.24385683D-07 
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5 
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5 

5 
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5 

5 

5 

5 

5 

5 

5 
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5 

5 

5 

5 

5 

5 

5 

5 


20 

0 

Yij  (a=  5) 
0.72711274D+00 

(CONT'D 

5  29 

.) 

4 

0.153014650-13 

20 

1 

0.25763815D+00 

5 

30 

0 

0.953790920+00 

20 

2 

0.15155124D-01 

5 

30 

1 

0.459954470-01 

20 

3 

0.939838530-04 

5 

30 

2 

0.213599460-03 

20 

4 

0. 609642090-08 

5 

30 

3 

0.351064440-07 

21 

0 

0.768079860+00 

5 

30 

4 

0.353110740-14 

21 

1 

0.221770630+00 

5 

31 

0 

0.961760910+00 

21 

2 

0.101056410-01 

5 

31 

1 

0.381020560-01 

21 

3 

0.438644870-04 

5 

31 

2 

0.137019180-03 

21 

4 

0.150617460-08 

5 

31 

3 

0.155471420-07 

22 

0 

0.803822080+00 

5 

31 

4 

0.810846880-15 

22 

1 

0.189461980+00 

5 

32 

0 

0.968396450+00 

22 

2 

0.669562570-02 

5 

32 

1 

0.315158410-01 

22 

3 

0.203175010-04 

5 

32 

2 

0.877015390-04 

22 

4 

0.368176010-09 

5 

32 

3 

0.686218760-08 

23 

0 

0.834720650+00 

5 

32 

4 

0.185336430-15 

23 

1 

0.160858010+00 

5 

33 

0 

0.973910730+00 

23 

2 

0.441199890-02 

5 

33 

1 

0.260332450-01 

23 

3 

0.934637130-05 

5 

33 

2 

0.560199130-04 

23 

4 

0.891373490-10 

5 

33 

3 

0.301936270-08 

24 

0 

0.861230940+00 

5 

33 

4 

0.421800150-16 

24 

1 

0.135871400+00 

5 

34 

0 

0.978485660+00 

24 

2 

0.289338880-02 

5 

34 

1 

0.214786270-01 

24 

3 

0.427270430-05 

5 

34 

2 

0.357144950-04 

24 

4 

0.213929640-10 

5 

34 

3 

0.132462370-08 

25 

0 

0.883834250+00 

5 

34 

4 

0.956080340-17 

25 

1 

0.114274280+00 

5 

35 

0 

0.982275620+00 

25 

2 

0.188952590-02 

5 

35 

1 

0.177016500-01 

25 

3 

0.194215680-05 

5 

35 

2 

0.227281960-04 

25 

4 

0.509356280-11 

5 

35 

3 

0.579522870-09 

26 

0 

0.903005880+00 

5 

35 

4 

0.215889110-17 

26 

1 

0.957638990-01 

5 

36 

0 

0.985411110+00 

26 

2 

0.122934450-02 

5 

36 

1 

0.145744550-01 

26 

3 

0.878197060-06 

5 

36 

2 

0.144394350-04 

26 

4 

0.120393300-11 

5 

36 

3 

0.252882710-09 

27 

0 

0.919195430+00 

5 

36 

4 

0.485750490-18 

27 

1 

0.800070330-01 

5 

37 

0 

0.988001970+00 

27 

2 

0.797136870-03 

5 

37 

1 

0.119888680-01 

27 

3 

0.395189710-06 

5 

37 

2 

0.915887840-05 

27 

4 

0.282662540-12 

5 

37 

3 

0.110078360-09 

28 

0 

0.932815820+00 

5 

37 

4 

0.108925870-18 

28 

1 

0.666686990-01 

5 

38 

0 

0.9901^0440+00 

28 

2 

0.515307980-03 

5 

38 

1 

0.985376250-02 

28 

3 

0.177045230-06 

5 

38 

2 

0.580068870-05 

28 

4 

0.659545920-13 

5 

38 

3 

0.478054570-10 

29 

0 

0.944238190+00 

5 

38 

4 

0.243481350-19 

29 

1 

0.554295360-01 

5 

39 

0 

0.991903690+00 

29 

2 

0.332192680-03 

5 

39 

1 

0.809264590-02 

29 

3 

0.789894690-07 

•5 

39 

2 

0.366857220-05 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 
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5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 
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39 

3 

0.20715698D-10 

5 

49 

2 

0.352051250-07 

39 

4 

0.54261559D-20 

5 

49 

3 

0.437909870-14 

40 

0 

0.99335618D+00 

5 

49 

4 

0.143130700-26 

40 

1 

0.66415059D-02 

5 

50 

0 

0.999107990+00 

40 

2 

0.23170054D-05 

5 

50 

1 

0.891985800-03 

40 

3 

0.89581398D-11 

5 

50 

2 

0.220032060-07 

40 

4 

0. 12058124D-20 

5 

50 

3 

0.186344630-14 

41 

0 

0.99455163D+00 

5 

50 

4 

0.311153700-27 

41 

1 

0.54469118D-02 

5 

51 

0 

0.999272120+00 

41 

2 

0. 14615010D-05 

5 

51 

1 

0.727868840-03 

41 

3 

0.33661445D-11 

5 

51 

2 

0.137407780-07 

41 

4 

0.26723410D-21 

5 

51 

3 

0.791964660-15 

42 

0 

0.99553471D+00 

5 

51 

4 

0.675269730-28 

42 

1 

0.44643646D-02 

5 

52 

0 

0.999406270+00 

42 

2 

0.92074797D-06 

5 

52 

1 

0.593717890-03 

42 

3 

0. 16654161D-11 

5 

52 

2 

0.857424630-08 

42 

4 

0.59072802D-22 

5 

52 

3 

0.336180920-15 

43 

0 

0.99634254D+00 

5 

52 

4 

0.146308440-28 

43 

1 

0.36568766D-02 

5 

53 

0 

0.999515880+00 

43 

2 

0.57939850D-06 

5 

53 

1 

0.484111720-03 

43 

3 

0.71612893D-12 

5 

53 

2 

0.534629490-08 

43 

4 

0. 13026310D-22 

5 

53 

3 

0.142540710-15 

44 

0 

0.99700588D+00 

5 

53 

4 

0.316503970-29 

44 

1 

0.29937564D-02 

5 

54 

0 

0.999605400+00 

44 

2 

0.36419377D-06 

5 

54 

1 

0. 394598780-03 

44 

3 

0.30741144D-12 

5 

54 

2 

0.333115310-08 

44 

4 

0.286578820-23 

5 

54 

3 

0. 603701830-16 

45 

0 

0. 997550200+00 

5 

54 

4 

0.683648590-30 

45 

1 

0.244957570-02 

5 

55 

0 

0.999678470+00 

45 

2 

0.228679990-06 

5 

55 

1 

0.321526210-03 
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